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PART  I:  ULTRASONIC  ATTENUATION  IN  NORMAL  METALS  AT  LOY  TEMPERATURES. 

Below  about  20°  K  the  attenuation  of  sound  waves  in  a  single 
metal  crystal  arises  primarily  from  the  electron-ion  interaction. 

’’lien  a  sound  wave  passes  through  a  metal  the  ions  are  displaced  from 
their  equilibrium  positions  and  an  electric  field  is  set  up  in  the 
metal,  Thi3  field  exerts  a  force  on  the  ions;  in  addition  there  is 
a  force  due  to  the  transfer  of  momentum  to  the  ions  in  electron-ion 
collisions.  Because  of  the  finite  relaxation  time  for  the  electrons, 
the  resultant  force  is  out  of  phase  with  ionic  motion  and  hence  leads 
to  attenuation  and  dispersion  of  sound  waves.  The  problem  in  formu¬ 
lated  in  terms  of  the  Boltzmann  transport  equation,  '"hen  the  colli¬ 
sion  integral  in  the  Boltzmann  equation  is  characterized  by  a  single 
relaxation  time  our  procedure  gives  the  same  results  as  obtained  by 
previous  authors. 

Next  a  method  for  solving  the  Boltzmann  equation  at  low  tem¬ 
peratures  is  given  when  the  collision  integral  is  due  to  the  interac¬ 
tion  of  the  electrons  with  phonons  and  impurities.  Although,  in 
principle,  the  method  is  applicable  for  all  relevant  wave-lengths, 

A  ,  of  the  sound  waves,  explicit  solutions  have  been  obtained  only 
for  the  limiting  cases  AXX?  and  A<A  4  ,  where  /  is  the  mean  free  path 


(in) 


of  the  electrons.  5'orA>>i?  it  is  found  that  the  relaxation  time, 

*  appropriate  for  attenuation  is  the  same  for  both  shear  and  di- 

* 

lational  waves.  However,  T#  is  different  from  Zh  y  the  relaxation 
time  appropriate  for  electrical  conductivity,  .  farther,  it  is 
found  that  the  relaxation  time  appearing  in  the  expression  for  the  at¬ 
tenuation  for  extremely  high  frequency  shear  waves  is  to  be  interpreted 
as  Zc  o  Lastly,  it  is  shown  that  where  electron-electron  collisions 
are  important  their  relaxation  time  may  be  deduced  from  the  observed 
temperature  variation  of  the  attenuation  and  cr, 

PART  II:  ULTRASONIC  ATT  MUTATION  IN  POLYCRYSTALS. 

In  a  polycrystal  the  crystallographic  axes  of  the  different 
grains  are  differently  oriented  with  respect  to  a  fixed  coordinate 
system.  Since  the  effective  elastic  constants  of  a  grain  depend  upon 
its  orientation,  they  differ  from  grain  to  grain  and  hence  the  sound 
waves  are  scattered.  In  this  part  of  the  thesis  we  calculate  tho  scat¬ 
tering  (and  hence  the  attenuation)  of  plane  3hear  and  dilational  waves 
due  to  the  random  orientation  of  the  grains  in  a  polycrystal.  These 
calculations  are  a  generalisation  of  the  work  of  earlier  authors  to 
take  into  account  the  fact  that  an  elastically  anisotropic  medium  has 
actually  several  independent  elastic  constants.  The  final  expressions 
for  the  attenuation  coefficients  exhibit  qualitatively  similar  features 
as  the  corresponding  expressions  of  previous  work  but  differ  quantita¬ 
tively  from  them. 
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A  plane  shear  or  a  plane  dilation-  1  sound  wave  propagating 
in  solids  suffers  at tenu;  tion.  T1  e  attenuation,  o(  ,  i 

reciprocal  of  the  distance  in  whicv  t  :e  intensity  of  the  w  ve  is  reduced 
to  l/e  of  its  value. 

In  the  case  of  nets Is,  whose  1 in err  dimensions  ore  large  com¬ 
pared  with  the  wavelength  A  of  t"  e  sound  ave,  t  e  tt(  nu  va¬ 

ries  roughly  v/ith  temperature  as  in  figure  1.  The  attenuation  at  room 
temperatures  and  higher  and  the  maximum  around  100°K  are  ue  to  the  in¬ 
teraction  of  the  sound  wave  with  dislocations.  T’  e  precise  temperature 
at  which  this  maximum  occurs  depends  on  t  gular  Freqi  -ncy .  ur  .  of 
t:ie  sound  wave  and  on  the  crystal.  The  theories  of  this  irt  -r&ction 
have  "been  developed  by  Koehler,  (1952)  [l],  and  Graneto  and  Tucke, 

(1956)  [2],  for  the  room  temperature  region  and  by  Seeger.  ^195 6  [ 3  ] , 

for  the  so-called  Ttordoni  peel:  at  intermediate  temperatures,  w-  shell 
not  be  concerned  with  these  portions  of  the  temperature  range  in  this 
thesis,  but  shall  deal  with  the  attenuation  below  about  2 O^TT.  Here 
again,  the  precise  temperature  dependence  of  the  attenuation  depend s 
upon  the  frequency  of  the  sound  wave  and  on  the  crystal. 

The  mechanism  producing  the  attenuation  belo  about  20°Y 
is  largely  the  interaction  of  the  sound  wave  with,  the  electrons,  as 
first  concluded  by  BSmrael  and  Mason,  (1956)  [4].  T^e  main  evidence 
for  this  conclusion  is  two— fold:  (l)  In  normal  net  Is  and  for  sound 
waves  with  wavelengths  much  greater  than  the  mean  free  path  of  th<  elec- 
tr  n  '  the  tei  ;  r-  ture  dependence  of  o(  Is  \  alj  tively  t)  e  same  s  ths  t 
for  the  el* ctrical  conductivity  .  T  is  shows  that  the  a  o  interaction  i 
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Dependence  of  Ultrasonic  Attenuation  in  Lead  on  Temperature, 

(!)  Per  the  Ncr/val  State  of  the  Heta-h  and  (2)  fa-  the  S  apef  couch  act  Strafe 
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to  a  large  extent,  responsible  for  both  of  these  effect  s  (2)  Tn  up  r- 
conducting  materials  the  attenuation  drops  sharply  at  the  transition 
temperature  to  a  very  low  value  and  decreases  rapidly  to  sere  as  the 
temperature  approaches  absolute  zero.  doth  shear  one!  dilation '  1  sound 
waves  behave  similarly  and  a  typical  plot  is  illustrated  in  ••Imre  2, 

Part  I  of  the  present  thesis  is  a  contribution  to  the  theore¬ 
tical  study  of  the  attenuation  of  sound  waves  in  normal  metals  result¬ 
ing  from  their  interactions  with  the  electrons  in  the  metal.  e  rela¬ 
tionship  of  our  work  to  that  of  previous  authors  will  bo  discussed  in 
the  introduction  to  Part  I. 

In  the  case  of  polycrystalline  materials  then  -'re  additio¬ 
nal  mechanisms  causing  attenuation.  One  of  these  arises  from,  the-  f  ct 
that  the  crystallographic  axes  of  the  different  grains  in  the  poly- 
crystal  ore  differently  oriented  with  respect  to  a  sot  of  coord In-  t 


axes  fixed  in  space,  that  is,  the  labor.'  tory  coordinate  system.  'a 
a  result  the  effective  elastic  constants  determining  the  propagation 
of  the  sound  wave  in  a  given  direction  varies  from  grain  to  grain  and 
this  gives  rise  to  o  scattering  of  the  sound  wave  .  When  ibt  linear 
dimensions  of  the  gr.  in  are  small  compared  to  the  wavelength  of  the 
sound  waves,  the  attenuation  arising  from  t  in  scattering  process  is 
given  by 


CV  -  Tif r 


(0 


where  7)  is  a  constant  depending  upon  the  specimen  used  end  it.r  elastic 
constants  and  f  *  *r/21T  is  the  frequency  of  tl  e  nd  wav  .  h<  va¬ 

riation  of  the  constant  P  v/ith  temperature  arises  from  the  variation  of 
the  elastic  constants  with  temperature  find  over  a  range  of  0JTC  to  yO^Y 


« 


Attenuation  Difference  in  Nepers  per  Centimeter. 
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Figure  2. 

Dependence  of  Ultrasonic  Attenuation  in 

Lead  on  Temperature  Below  IZ°K. 

O')  For  the  Mcr/na!  State  of  the  Metal  ,  (2)  /<>/*  the.  Sitter  - 

Conduct t  hj  State  .  /he  thermal  State  to  fiehounech  Belcoo  The 

Transit,  ion  Temperature  ,  Aj  }  /y  The  ft pph feat  ion  of  a  /toffietic  hje/oh- 
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is  not  more  than  JO;-,  In  view  of  the  frets  that  the  other  meet,  nistns, 
some  of  them  being  t  i.e  same  as  those  responsible  for  attenuation  in 

St 

sin--  le  crystals,  depend  sensitively  on  temperature  and  ’  ;  ve  a  frequency 
dependence  which  is  sub stanti ally  different  from  a  fourth  power  law,  it 
is  not  difficult  to  isolate  the  attenuation  arising  from  the  scattering 
due  to  random  orientations  of  the  grains,  ~a.rt  II  of  this  thesis  is 
a  calculation  of  I)  in  (i).  Again  the  relationship  of  our  work  i+h  that 
of  previous  authors  will  be  discussed  in  the  introduction  to  fart  IT, 

It  would  be  useful  to  note  the  frequency  ranges  and  wavelengths 
of  the  sound  waves  for  .hich  the  mo  chan  isms  studied  ''ere  contribute  to 
attenuation,  for  TVrt  I  the  first  measurements  were  made  in  the  mega¬ 
cycle  region,  but  very  recently  Lax,  (1959)  [5]»  has  extended  the  lower 
limit  to  close  to  ten  kilocycles  per  second.  The  upper  limit  to  the  ob¬ 
servable  frequency  range  is  determined  by  the  experiment  1  difficulties 
in  obtaining  very  high  frequency  sound  waves  of  appreciable  energies. 

In  Part  II  the  fourth  power  dependence  becomes  meaourcble  in  the  mega¬ 
cycle  range.  Taking  on  average  velocity  of  3  x  10^  centimeters  nor  se- 
cond  for  the  sound  wave,  frequencies  of  10  cycler  per  second  represent 

a 

wave  1  eng t'  s  of  about  30  centimeters  and  frequencies  of  10  cycles  per 
second  represent  wavelengths  of  about  0,003  centimeters, 

-ince  the  experimental  procedure  is  outside  th;>  scope  of  this 
thesis,  it  will  be  mentioned  only  Ly,  h  u  - 

rating  ultrasonic  waves  in  solids  is  that  of  piezoelectric  and  magnet o- 
ntrictive  transducer  - ,  such  as  quarts.  In  the  megacycle  r  n;  o  the  vulse 
technique  is  used.  This  technique  consists  of  sending  ultrasonic  pulses 
short  with  respect  to  the  length  of  the  specimen  but  of  about  one  hundred 
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wavelengths  no  as  to  ensure  ste-  dy  state  conditions  t1  rough out  no r t 
of  tho  pulse,  through  the  sample  at  regular  intervals.  the  velocity  and 
attenuation  of  the  initial  pulse  or  the  reflected  pulse,  (echo-t  jchni— 
sue),  are  measured.  A  block  diagram  of  the  experimental  apparatus  is 
shown  in  figure  3*  At  lower  frequencies  the  pulse  toe  ni quo  cannot, 


of  course,  he  used,  so  that  here  the  velocity  and  attenuation  are  mea¬ 
sured  hy  setting  up  standing  waves  in  the  sample  and  monnuriiv  their 


time  rate  of  decay 


Pulse 


7 


o 

< 


CO 

O 

U 

z 

bD 

•H 


O 

■P 

© 

m 

D 

© 

•p 

aj 

U 

cn 

Cu 

£X 


© 

4-^ 

c 

© 

p 

*—4 

•H 

M 

© 

Pu 

K 

W 

© 

£ 

4-> 

O 

P 
*•— • 

aj 

P 

bO 

aJ 

•H 

o 

o 

o 

I — I 

PQ 


Study  Sound  Waves  in  Solids. 


PART  I 


ULTRASONIC  ATTENUATION  AND  DISPERSION 
IN  NORMAL  METALS 
AT  LOW  TEMFETIATURES . 


8 


1 


T^TrrrM  /A-\TT^xrn  t  fYM 

•  .l  i  o  X  x. 

As  air  >ad;  mentioned,  the  attenuation  of  round  ?/. -vos  in  metals 
at  liquid  helium  temperatures  arisen  from  the  interaction  "between  the 
conduction  electrons  and  the  ions  in  the  metal.  In  this  section  wo  shall 
review  "briefly  the  earlier  theoretical  work  and  then  outline  the  work 
done  in  this  thesis. 


1.1  hBfUI.tfi  Crl  hh  lf  "L  W-TdK 

The  first  t'hocrc tier-1  work  on  the  subject  is  duo  to^Mason  and 
iSrnrnel ,  (195  ^  J 4 3»  w^°  ascribed  the  attenuation  to  bo  due  to  the  vis¬ 
cosity  of  t-ie  electron  gas  in  the  metal.  If  the  electron  gas  is  asrumed 
to  be  rigidly  coupled  to  the  ions,  the  viscous  force  per  unit  area  on 
an  clement  of  volume  of  the  metal  contributes  to  the  shearing  stress 
a  term  of  the  type 


•\r 


T 


/Z 


yU  5, 


a 


+ 


7 


(i) 


♦ 

here  S,z  is  the  shearing  strain,  siZ  is  its  time  rate  ..of  cV  nge%  y  is 
the  coefficient  of  viscosity,  and  u  is  the  1  e; r  elastic  constant.  Th« 
attenuation  coefficient,  oC  ,  is  then  obtained  in  t’  e  u  ual  way  ( Seo 
for  example,  !*ason ,  loctromechanical  Transducers  anu  v/evo  Alters) , 
resulting  in 


c*r 


s 


(2) 


for  shear  waves  and 


cx 


^s(fr 


X 


(?) 
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for  dilational  waves, 
the  velocities  of  the 
tivoly,  in  the  metal, 
types  of  waves,  has on 
for  the  electron  gas. 


X  is  the  eompressional  viscosity,  V  and  v.  are 
shear  waves  ind  of  the  dilational  waves,  respoc— 
from  the  experimental  evidence  or.  o(  for  the  two 
and  BSmmel ,  J*4l  »  haye  coj  t  t  —  < -O. 


? 


The  kinetic  theory  expression  for  the  viscosity  of  a  gas  is 
=  1  Nra/v,  (4) 

3 


where  I»  is  the  number  of  particles 
particle,  /  is  the  mean  free  path, 
trons  we  have 


per  unit  volume,  m  is  the  mass  of  a 
and  v  is  the  mean  velocity.  For  elec 


ltnv*  = 

2  10  m 

where  II  is  now  the  number  of  free  electrons 
that  expressions  (2)  and  (f)  for  attenuation 
square  of  the  frequency.  The  temperature  do¬ 


per  unit  volume.  We  note 
are  proportional  to  the 
vend once  of  the  attenua¬ 


tion  comes  about  through  the  dependence  of  J  on  temper  turef,  •■Jason 
and  BSmmel  assumed  that  the  mean  free  path  appearing  in  (4)  may  be 
taken  to  be  the  same  as  that  which-  appears  in  the  expression  for  tho 


electrical  conductivity,  cr ,  when  written  in  the  form 


cr* 


ITe*  f 
m  v 


or  J?  -  c r  mv 


(0 


.  Q' 


T) us  the  attenuation  would  be  proportional  to  the  electrical  conducti¬ 
vity. 

A  treatment  such  as  above,  in  terns  of  the  viscosity  of  the 
electron  gas,  is  valid  only  as  long  as  the  mean  free  path  is 


small  com- 
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pored  to  tho  wavelength,  _/\.  ,  of  the  sound  wave ,  that  is.  for  kJ?  <  1 
where  k  =  2  TT /A  »  W  en  />A  it  iq.  no  longer  permi  ible  to  i  ;ol,  to  for 
consideration  a  portion  of  the  metal  small  compared  to  A  ,  Moreover, 
the  passage  of  a  round  wave  through  the  metal  producer,  an  electric  field 
in  the  metal  and,  therefore,  in  the  presence  of  tl  e  sound  wave  the  elec¬ 
trons  are  affected  both  by  collisions  with  the  ions  and  by  the  electric 
field.  A  comprehensive  treatment  which  takes  hot-  these  effects  into 
account  lias  been  given  by  Pippard,  (1955)  [  6 .1 ,  in  terms  of  the  free  elec¬ 
tron  model  of  a  metal. 


The  basic  ideas  underlying  his  work  are  a:  follows.  When  the 
lattice  ions  .are  displaced  from  their  equilibrium  po  i lions  due  to  an 
ultrasonic  -save  passing  through  a  metal,  an  electric  field  in  set  up 
inside  the  metal.  This  electric  field  and  the  rtrtionary  di  tribution 
function  of  the  electrons  are  determined  by  the  collision.;  of  the  elec¬ 
trons  with  the  lattice  vibrations  and  other  defects  in  the  metal  ,  by  the 
equation  of  continuity,  and  by  I.taxwell's  electromagnetic  equations,  "ho 
resultant  distribution  function  is  not  identical  with  that  whicl  the 
electrons  would  have  if  they  were  locally  in  thermal  equilibrium.  Hence 
there  is  a  dissipation  of  sound  energy  into  heat  energy.  T'ippard  adopted 
the  kinetic  approach  to  the  problem,  namely,  that  of  following  n  single 
electron  at  the  formi  surface  and  considering  the  changes  in  its  velocity 
due  to  both  the  electric  field  generated  by  the  sound  w  vea  and  the  col¬ 
lisions.  for  the  latter  he  assumed  that  th  velocity  rela  :os»  exponen¬ 
tially  to  its  value  corresponding  to  local  equilibrium  with  a  relaxa- 
n  time,  Z  •  The  attenuation  n  ultr  nic  i  tl  n  obtained 

by  calcul' ting  the  Joule  heat  associated  with  the  .  t -a  :y  state  electronic 


current  act  up  by  the  wave 


Tiring  the  same  be  .sic  ideas  as  those  used  by  pippord  St.  ;'nborg 
(1958a)  [7],  has  formulated  the  problem  in  terms  of  the  Maxwell—' Boltz- 
mann  transport  equation.  In  this  work  the  collision  int  gr  1  is  e  mpli- 
fiod  by  the  assumption  that,  due  to  the  collisions,  the  di  tribution 
function,  f,  for  the  electrons  relaxes  towards  a  distribution  function 
f,  which  the  electrons  would  have  if  they  were  locally  in  thermal  equi¬ 
librium,  with  a  relaxation  time,  Z  ,  so  that 


|_at 


coll. 


(f  -  h 

L 


7) 


f  t  rinberg  then  sets  up  the  equations  of  motion  for  the  ions  by  adding 


to  the  elastic  stress  tensor  a  kinetic  stress  tensor  arsocirted  v.  it} 


n 


the  st 9  dy  state  distribution  function  of • the  electrons.  The  components 
of  this  kinetic  stress  tensor  are  given  by 


vy 


(2 s| j/  Vy  [fd0  -?(»)] 


dv 


'0) 


'"bis  procedure  is  somewhat  analogous  to  th?  t  used  in  the  kinetic  theory 
of  gases  in  order  to  obtain  the  Tavier-Stokes  equation  of  hycirodym  mice 
from  the  I.Iaxw el 1- Bol t zmann  transport  equation.  By  solving  Ms  onus- 
tionn  of  motion  for  plana  waves,  Steinberg  is  able  to  calculate  the  de¬ 


pend  enco  of  both  the  attenuation  and 


the  sound  velocity 


on  frequency.  The  analysis  is  carried  out  for  both  sheer  -nc  dilotio— 
nn.l  waves  snd  the  results  on  attenuation  are  in  agreement  v  -t.b  tho  o 
obtained  by  "ippard.  However,  teinb  rg’s  own  conclusion  that  his-  re¬ 
sults  are  at  variance  with  T’ippard’s  is,  in  our  opinion,  based  on  n  er¬ 
roneous  argument  am  we  discuss  this  in  section  2.4  of  this  thesis. 
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Independently  of  T  ippard  and  3t*  inberg,  ‘kheizer  ct  £l  , 

(1957)  [0],  have  also  made  a  substantially  similar  Boltznr  nn  formulation 
of  the  problem.  Their  analysis  consists  of  •  calculation  of  the  at  term- 
ation  for  dilational  waves  in  the  limits  of  very  low  an  very  high  fre¬ 
quencies  by  calculating  the  entropy  change  in  the  electron  gas. 

In  the  "forks  described  above  the  various  authors  have  assumed 
tl  e  existence  of  s.  relaxation  time  for  the  distribution  function  of  the 
electrons.  [Oven  if  a  relaxation  time  for  the  attenuation  problem  does 
exist,  it  is  not  evident  that  this  relaxation  time  should  bo  the  some  as 
that  appearing  in  the  expression  for  electrical  conductivity  r.s  is  assumed 
BSmmel*  Nor,  as  pointed  out  by  Steinberg  ,  (1958?  ) >  is  ther 


any  "a  priori"  reason  to  assume  that  the  relaxation  times  appearing  in 
the  expressions  for  the  attenuation  of  dilation? 1  and  shear  waves  are 
the  same.  There  is  also  the  question  as  to  whether  tl  o  same  relaxation 
time  is  appropriate  to  describe  the  attenuation  at  both  very  low  and 


rery  high  frequencies  in  the  cs.se  of  shear  waves. 

Before  outlining  the  work  done  in  the  present  thesis,  we  note 
that  -te  inberg,  \  1958b)  j_9j»  has  calculated  the  viscosity  of  a  degene¬ 
rate  electron  gas.  lie  finds  that  at  low  temperatures  the  relaxation 
time  resulting  from  the  electron-phonon  interaction,  appropriate  for 
e  viscosity,  is  approximately  one-third  the  relaxation  time  appro¬ 
priate  for  electrical  conductivity.  Wo  shall  designate  this  latter 
relaxation  time  by  2T,  hereafter.  Tor  impurity  scattering  one  obtains 
a  relaxation  time  £  times  Z*/  ,  the  residual  relaxation  time  "or  elec¬ 
trical  conductivity,  where  £  is  a  constant  of  the  order  of,  but  lose,  than, 
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unity.  These  results  are  relevant  for  the  attenuation  of  lev; 

cy,  ki,<^l,wzr/^<l,  shear  waves. 


•u 


froquen 


1.2  OUTLIE-:  T". '  WOPP  IV  TPI.  "'”E  T  . 

In  the  present  work  we  carry  out  several  progressive  c*  Icu- 
lations.  These  calculations,  in  common  with  those  described  previously 
use  the  hommerfeld  free  electron  model  of  a  metal,  '  yen  or  lization  to 
&  model  where  the  electrons  are  do.  cribed  by  a  sing  le  effective  ness 
throughout  the  band  is  straightforward.  We  also  use  in  our  ca.lcul  - 
tions  Pippard’s  basic  idee  that  the  electric  field  produced  by  the 
passage  of  a  sound  wave  and  the  stationary  distribution  function  of 
the  electron  are  determined  by  the  various  kind3  of  collisions  the 
electrons  can  suffer,  by  the  equation  of  continuity,  and  by  Maxwell1 s 
electromagnetic  equations » 

In  section  2,1  we  formulate  the  problem  in  terms  of  the 
Yaxwe  11— Boltzmann  transport  equation.  The  first  step  is  to  as-  ume  that 
the  collision  integral  can  be  characterised  by  a  unique  relaxation  time, 
r  ,  in  accordance  with  (?)•  We  are  now  able  to  determine  th  electric 
field  and  the  stationary  distribution  function,  f ,  for  the  electrons  in 

terms  of  T .  li  art  of  the  work  is,  of  course,  identic;  1  with  Ptein— 

* 

berg's  procedure.  Wo  next  set  up  the  equations  of  motion  for  wave 
propagation  in  a  metal.  Instead  of  calcul  ting  tv'e  kinetic  stress  tensor 
t  f  as  Steinberg  has  done,  re  use  the  folio  in 
ideas.  The  electric  field  set  up  in  the  metal  due  to  the  pas  age  of 

■?  ' 

Our  work  v/oo  done  irndg.  ndentlj  of  Steinberg's  and  summary  was  sent 
to  tie  "International  Conference  on  the  Transport  Properties  of  -vtals 
at  Low  Temperatures”,  Geneva,  Pew  York,  1P5^«  ome  of  the  results  on 
the  calculation  of  the  relaxation  times  contained  in  t'  is  p.  rt  of  this, 
thesis  were  reported  at  the  January,  I960,  meeting  of  the  Am  ricon  Phy¬ 
sical  Society,  Bulletin  of  the  American  Physical  Society  i960  Annual 
Meeting  at  few  York,  C12,  Page  14. 
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the  ultrasonic  wave  exerts  a  force  on  the  positively  charged  iorr ,  the 
electric  force.  rbrther,  in  electron— ion  collision."  there  if?  »  direct 
transfer  of  momentum  from  the  electrons  to  the  ions  per  unit  time  end 
this  results  in  a  force  on  the  ions,  the  impulsive  ^orce.  Flnrlly  there 


is  the  elastic  restoring  force  between  the  ions  indepench-nt  o"  the  elec¬ 


tron-ion  interaction,  that  is,  even  when  the  impul  iive  1  t)  1  c— 

trie  force  are  zero.  The  resultant  of  these  three  force''  d  ,  because 
of  the  finite  relaxation  time  for  the  electrons,  in  general,  net  in 
fhase  with  the  displacements  of  the  ions  and  hence  the  waves  are  atte¬ 
nuated.  In  this  manner  we  are  able  to  calculate  both  the  a ttonuation 
and  the  dispersion  of  the  sound  waves. 

The  expressions  for  the  attenuation  are  identical  with  those 
obtained  by  Pippsrd  and  teinb  ?ry.  However ,  our  procedure  of  setting 
up  the  equations  of  motion  for  the  metal  gives  us  an  interesting  insight 
into  the  significance  of  z  in  one  particular  case.  Tt  i r-  'oond  that  the 
attenuation  of  a  shear  wave  with  an  angular  frequency  such  that  *rz>y  1 
results  solely  from  the  impulsive  force.  Since  this  impulsive  force 
arises  from  the  direct  transfer  of  momentum  from  the  electrons  J  t’  e 
ions,  it  will  contain  the  same  relaxation  time  as  that  appearing  in 
the  expression  for  the  electrical  conductivity.  hence,  in  this  high 
.frequency  limit,  the  attenuation  dll  also  contain  this  relaxation 
time.  A  confirmation  of  this  expectation  is  obtained  in  section  3.5. 

che*j, es  Zn  fde.  t/e/ocity  y/  t/tc  sound  tjaoes 

As  regards  to  the- dispersion,  : •  na  1 1  o o i ty  eh »  e y >■] ,  of 

the  order  of  m/m ,  the  mass  of  the  electron  over  the  mass  of  t’  ion, 
are  associated  with  the  attenuation.  Incidentally ,  our  methods  also 
yield  an  expression  for  the  contribution  to  the  velocity  of 


the  di.li 
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tional  waves  arising  from  the  v.  riation  in  the  electronic  chrrg  den¬ 
sity  associated  with  the  passage  of  this  type  of  wave  in  e  mats 1 . 

This  contribution  is  frequency  nt  f  r  high  fre  •  enc  oun 

to  v  »* 

(  >  10  radians  per  second),  but  is  independent  of  the  relaxation 

time  for  the  electrons.  e  discuss  this  contribution  furt’  er  in  sec¬ 
tion  2.4. 

In  sections  3  and  4  we  discard  the  assumption  of  t:  e  oxi.  tence 
of  a  single  relaxation  time  and  consider  the  solutions  o’  td  ■  transport 
equation  when  the  collision  integral  is  due  to  the  elect  -on-phonon  in¬ 
teraction,  (section  3 ; »  and  to  the  electron-impurity  scatte  ring  ^sec¬ 
tion  4).  ?or  the  electron-phonon  interaction  we  start  itr  th<:-  col¬ 
lision  integral  as  given  by  Wilson.  First  ve  modify  it  suitably  to 
take  into  account  the  fact  that,  in  the  presence  of  a  sound  v? v-  tv 
distribution  function  for  the  electrons,  f,  corresponding  to  loo  1  equi¬ 
librium  is  different  "rom  the  equilibrium  distribution  Emotion  f 

^  o ' 

in  the  absence  of  the  sound  wave,  f  is  just  the  usual  ferri-Di  me  '.’is- 

•  ^ 

tribution  function •  Secondly,  in  order  to  obt  in  I 

Foltzmann  transport  equation  we  expand  f  in  a  series  of  spheric  1  har¬ 
monics, 

f  «  f  +  HVc'’  (v)Y"  3fu  .  r9\ 

go  r  r  aTr 

ffe  find  that  it  is  necess  ry  to  retain  only  the  n=0  ter’  for  incident 
in e  dilations  1  waves  and  the  n=l  terms  for  incidenl  plan<  sh  ar  rsves, 
TU-ina  (9)  in  the  transport  equation  results  in  an  infinite  s^t  on  cou¬ 
pled  linear  integral  equations  for  the  c”(v) ’  s.  Tender  certain  simply- 
fying  assumptions  vs  lid  at  low  temperatures,  it  is  found  t  -sob 
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C^(v)  relaxes  towards  its  value  corresponding  to  local  equilibrium 
with  a  relaxation  time,  2Tp  ,  such  that 


7* 


PIP  + 


IT 


'  > 


'10) 


where  Z(  is  the  re]  tion  time  for  C.(v)  L  just  that  appropri ste 
for  electrical  conductivity.  The  set  of  coupled  linear  intogr  1  cqu;  - 
tions  have  now  been  reduced  to  a  set  of  coupled  liner r  a.lgebr<: ic  equa¬ 


tions. 


The  infinite  set  of  coupled  linear  algebraic  equations  for 


n 


the  Cy(v)fs  can  be  readily  solved  for  certain  speci  1  c 
corresponding  expressions  for  the  attenuation  derived.  firsts  v?h  n  Ic^r 
«  k£  <  <  1,  one  finds  that  the  relaxation  time  appear in  in  tl  x]  •- 
cions  for  is  the  same  for  both  dilation&l  and  she.'  r  wav  ~:3 . 
found  that  this  relaxation  time  is  just  Z  »  T(  /3 .  Secondly,  for 
/,  >  1  but  wZt  <  1 ,  the  expressions  for  <*  ro’o  in  •  re¬ 
laxation  times,  fhirdly,  for  «rZ,  >^1,  the  expression  for  oi  of  i- 

lational  waves  is  just  that  of  the  previous  co.se  but  the  ex  free  ion 
for  the  attenuation  of  the  shear  wr.ves  contains  a  relaxation  time  s’ ic> 
is  now  to  be  interpreted  as  2T,  . 

For  the  intermediate  ca»r,e,  ^  1,  t'  e  solutions  of  the  cou¬ 
pled  set  of  equations  are  much  more  difficult  to  obtain  and  no  h;  ve 
not  found  s  way  to  obtain  analytical  solutions.  An  attempt  to  obtain 
numerical  solutions  of  these  equations  has  not  been  made  in  this  thesis. 

In  section  4  wo  discuss  briefly  the  case  when  the  collision 
integral  arises  from  impurity  scattering  alone  and  find  results  which 
are  qualitatively  the  same  as  those  obtained  in  section  ?  with  respect 
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to  relaxation  times.  In  particular  one  can  write  =*  t  ■  ■■  ^ 

depends  on  the  interaction  potential  between  the  electrons  and  tvs  im¬ 
purity  ions.  In  the  cane  of  3  simple  screened  coulomb  -potential,  e.p 
depends  only  on  the  range  of  the  interaction  poteiitial.  The  value  of 
&2  has  been  calculat'd  as  a  function  of  the  screening  distance  end  is 
found  to  lie  between  one-third  and  unity. 

Section  5  consists  o"  a  discussion  of  t’  e  r<  suits  obtr inod 
in  sections  2  to  4  in  relation  to  experiment.  The  influence  of  elec¬ 
tron-electron  collisions  is  brought  into  this  discussion  on  a  quali¬ 


tative  level 


2.  CALCULATION  0*  7 II I]  ATTENUATION  Ain  THE  DIOPHRUION  RELATIONS  A°- 


SUITING  THE  EXISTENCE  O'  A  OIHGLB  RELAXATION  '  l!'E 
2.1  FORMULATION. 

The  first  step  in  our  calculations  is  to  formulate  the  des¬ 
cription  of  the  electrons  in  a  metal  in  terms  of  the  "? xwel 1-Bolt zmnnn 
transport  equation. 

Let  f  (v)  denote  the  Fermi-Dirac  equilibrium  distribution 


function  for  the  electrons  in  the  undisturbed  metal,  where 

f„Cv)  '  '  1-1 

2 

2 

mvn  m  is  the  mass  of  the  electron,  v  is  its  velocity. 


[exp  [  ^  nv2  -  Eo  )/  kT}  +  l]  . 


(11) 


In  (11),  E 


o 


v^  is  the  velocity  of  an  electron  at  the  surface  of  the  Fermi  sphere, 
k  is  the  Boltzmann  constant,  and  T  is  the  t  -niperature.  '^he  number  of 
free  electrons,  ,  per  unit  volume  and  vq  are  related  by 


N 


/  2m 

(T73T 


i  r  tjDi? 


3 


(l?) 


It  becomes  convenient  to  ta2  e  the  number  of  free  electrons  to  be  equal 
to  the  number  of  motal  atoms  and  this  is  done  throughout. 

Let  an  ultrasonic  wave  of  wavelength,  A  ,  in  the  metal  be 
characterized  by  a  displacement  s  and  a  particle  velocity  if.  For  a 


plane  wave  travelling  in  the  x-direction. 


fexp  -  kx)]  ,  u 


* 

— * 

s 


(13) 


w) ere  k  »  2 Tf/Y.  The  distribution  function,  f,  for  the  electrons  in 


local  equilibrium  is  given  by 

f  = 


(oxp  [  (I  ‘  v*)ATj  +  • 


•'14) 


That  is,  T  is  a  -’Brrni  sphere  moving  with  the  lattice  and  having  a  surf  ce 
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1 


energy  Eq  +  The  additional  energy,  is  necessary  to  take  into 

account  the  density  fluctuations  of  the  electrons  caused  by  th  passage 
of  the  ultrasonic  wave  through  the  metal.  Here  we  assume  that  the  ef¬ 
fects  of  heat  flow  due  to  local  changes  in  temperature  of  the  electrons 
and  of  the  lattice  are  negligible.  e  note  that  the  attenuation  per 
wavelength  resulting  from  heat  flow  is  at  the  most  of  the  order 

Av  “  »  w^ere  Cp  and  Qv  are  the  specific  heats  at  constant  pres¬ 

sure  and  volume,  respectively.  At  liquid  helium  temperatures  this  fac- 

—8 

tor  is  of  the  order  of  10  or  less.  To  first  order  in  u  and  ,  '14) 
may  be  written  as 

f  *  fo  -  (m7,“  +  h  >  |4° '  (15 

o  Hi 

where  E  =  mv1' 

2 


Under  the  influence  of  the  ultrasonic  wave  a  certain  electric 

— * 

field,  £  ,  will  be  set  up  and  the  stationary  distribution  function  f. 
for  the  electrons  will  be  determined  by  the  Maxwell-Bo It zmann  transport 
equation 


+  v,  V7  f  +  ej  .  \j  f 
at  r  m  ^ 


where  e  is  the  electronic  charge. 

Bor  the  first  part  we  assume  that  f  relaxes  towards  f 
oue  to  collisions  with  a  relaxation  time,  t  ,  "  that  is, 


(H) 


coll 


_  I  .  (f  -  f) 
r 


Let  the  solution  of  (16)  be  of  the  form 

f  =  f  -  f  (v)  .  u  bf  , 

o  I  v  '  -~0  9 

air 


P°9 


(n) 


f  See  nex t 


e  1 for  -footnotes* 


(18) 


oo 


Substituting  (17)  and  ( 18 )  into  (l6)  and  retaining-;  only  term:'  of  first 
order  in  f,  3,  and  we  have 

1  (ezv.£  +  mv.u)  +  ^  ^  fl9^ 


fx  (v)  .  u 


TTmT?) 


(l  -  ia'cos s') 


whore 

+  a' 


-1 


kvr(l  +  i  ur  c')  , 


(?0) 


and  6*  is  the  angle  between  v  and  the  direction  of  propagation  of  the 
sound  wave,  (the  x-axis.  ). 


4f 

In  general,  there  will  be  also  an  induced  magnetic  field  produced 
but  in  the  first  order  theory  in  u  developed  here  this  will  have  a 
zero  effect  as  is  demonstrated  as  follows. 

This  term  will  be  of  the  form. 


e,(v  x  Tl)  .  \/v  f 
m 


e(v  x1)yu  \  fQ 

m  '• 

e(v  x  n).\7  f  + 
m 

/  ^  v  <4 

£(v  x  I!)  .mvjf c  * 
m  £  K 


-  t,  (v).  u>f  f 

as  ) 

higher  order  terms  in  u. 
0  , 


since  (v  x  r’W*  «=  o  and  since  T!  will  be  proportional  to  u.  Fence  the 
induced  magnetic  field  may  be  neglected  throughout. 


T  olstein,  (1959)  [l0l ,  has  given  a 
the  form  of  T  and  hi 3  treatment  will 

5.1. 


theoretical  justification  of 
be  discussed  more  fully  in  section 


The  electronic  current  density,  J  ,  and  the  number  of  elec- 
trons,  N,  per  unit  volume  at  a  point  in  space  are  given  by; 


N 


2m 

y 

h-' 


3 


f(v)dv*  and 


petrT  |  vf(v)dv 


(21) 


’he  equation  of  continuity  demands  that 


d  (N  ) 
«)  t 


div  J 


e 


0 


(22a) 


On  integrating  (14)  (throughout)  over  all  velocities  of  the  electron,  we 
see  that  (22a)  is  equivalent  to 


f (f  -  f)  dv 


0 


(  ’?b) 


The  electric  field,  £,  is  related  to  the  total  current  den- 

nity,  J  =  J  -  N  eu,  by  Maxwell’s  equations  (m  Gaussian  units), 

e  o 


curl  curl  t 


grad  div  £  - 

—  - 


N  eu)  + 
o  ' 


?]■ 


(  25) 


Equations  (22)  and  (23)  used  in  conjunction  with  fl9)  and 
(21)  are  sufficient  to  determine  ,  £  ,  and  £,  completely y  in  terms 

of  the  particle  velocity  u  associated  with  the  ultrasonic  wave. 

Next,  we  set  up  the  equations  of  motion  of  the  ions.  Ihe 
equations  of  motion  for  wave  propagation  in  the  x-direction  in  on  iso¬ 
tropic  solid  may  he  written  as 


/> 


V  e  e 

o 


-  I 


(°4) 


where  N^M)  is  the  ionic  density  of  the  metal,  -e  is  the  ionic 

chcrge,  and  a  is  the  square  of  the  velocity  of  the  wave  in  th  absence 

of  the  electron-ion  interaction,  that  is,  for  £  and  7  zero  in  (?d). 

The  second  term  on  the  right  hand  side  of  (24)  is  due  to  the  electric 

field  and  is  the  force  on  a  unit  volume,  the  electric  force.  The 

third  term  on  the  right  hand  side  of  (24)  is  the  impulsive  force  re- 

ulting  from  the  direct  collisions  of  the  ions  with  the  elections. 

The  term  I  can  be  determined  as  follows.  Due  to  the  col- 
c 

lisions  between  the  electrons  and  the  ions,  the  momentum  of  th'-'  elec¬ 
trons  per  unit  volume  changes  with  time  at  a  rate 


2m- 


J 


r  r 
mv 


G>f 

Ut 


dv 


coll . 


substituting  (17)  and  (21)  into  the  above  vie  v 


uve 


+ 


N  mu 
0 


(25) 


The  momentum  imparted  to  the  ions  per  unit  time  p  r  unit  volume  by  the 

electrons  in  collisions  with  the  ions  is  just  -  I  . 

c 

On  substituting  (13)  into  (24)  and  putting  k  =  -  ik0 

one  can  obtain  the  velocity,  V  *  ^/k^,  and  the  attenuation. o<  =  2k? , 
in  nepers/cm.  When  k2<<k^,  as  is  the  case  in  the  problem  considered 
here,  we  rewrite  (24)  in  the  form; 

?  p 

“  -  +  z 


2ik^kp )  + 


where 


nt\ 


-  (  Kef  +  7  )  S 

o _ c 

2 

yo  s 


Equating  real  and  imaginary  parts  in  (26)  then 


f'y 

c. 


7) 


oC 


and 

V2 


(’8) 


(29) 


2.2  INCIDENT  PLANK  PI  NATIONAL  WAVE::. 

We  represent  a  dilational  wave  propagating  in  th  x-dir action 
by  (13)  where  new  =  si.  3oth  u  and  £  will  be  along  the  x-cirec- 
tion  and  curl£  will  be  aero.  In  this  case,  (19)  h  conies 


fx  (v)  u 


(erv*^  +  rav^u)  +  $  t 

(l  4*  i*/^)  (l  -  ia'con^) 


(30) 


etting 

_  o 

dE 


-  £ ( E  -  Eq),  (21)  becomes 


e 


3NQeg 

2 

(l  4-  Wr)  mv 


|vq  (eT£  +  mu)  4  ia  ^ 


(31) 


and 


N  =  N  +  i 
o 


3NC^ 


2(1  +  iarr) 


•JL  i_i 


o 


vQ(e  4<f  4-  mu'  4  ia ^ 

1 

}  H-1 

g 

where 


£ 


1-1 


dx  and  a  = 


v  k  r 

o 


2  ^  1  -  iax 


1  +  i  ur  c. 

-  / 

and  x  ®  cos  6?  . 

Using  (22a)  we  may  now  solve  for  g  and  obtain 

£  =  i  v  g  (e  Z£  +  rru) 

->*  o°  _ *  * 

a(g  +  i^^T ) 

Using  (3l)  and  (34)  in  (23)  we  obtain  for  <£*  , 


(33) 


(34) 


e  £ 


. .. JL*.  r- JLL^ 
(1  -  i/X^/r) 


m 

4/TN  e2 
o 


26 


(32) 


where  ^ 


and 


"7 


x  - 


ii 


+  1^ 


r)  a 


1  + 


3s 


s  \ 

litre) 


k/a 

3g 


g  \.  00 


i  urc. 


Substituting  (31),  (34),  ana  (35 )✓  into  (25),  then 


I  -  IT  mu  i  Sur  (  X  -  1 ) 

C  0 


(i  +  iKsX/r) 


37) 


Hence,  from  (27) 


(K  e 
o 


>v  + 


I  )s 

c ' 


mN 


A  c 


(  X  -  1)  (1  +  i  ^/V) 
(  1  +  i  gf'^'X/Z") 


f  38) 


0 

Hirst,  v/e  note  that  #  appearing  above  i<?  the  reciprocal 
of  the  square  of  the  plasma  frequency  and  so 

2?  10  '  in  metals.  At  the  ultrasonic  frequencies  currently 
available  and  for  the  values  of  the  relaxation  time,  d  ,  that  occur 
in  metals,  the  terms  {/fa-jz  )  and  ( / or/-T  )  may  be  neglected,,  '"his 
is  demonstrated  in  the  following. 


fet 


A  “  X  y  +  i  %  2 


(39) 


Then  (38)  may  be  written, 

rf(*i  ■ x)  +  '( +  foi1 

V  -T-  L  r-X 


XP 


z2 


r 

%iy>  -  d"  *  K))x  [  1  -  2t'X+  x2+  *7 )] • 


(4o; 


Trom  (33), 


8 


- 1  r _ dx_ 

2  J  (  l-iaa 


and  therefore, 


+  1 


In 


2i 


<r> 


m) 


1  - 


_  i 

ltan  e 
a 


X  *  j*  _ a 

3  La-  tan 


a 


14/ ^  _ 


(4?) 


tan^  a 


For  a  <  <  1 ,  k,  / < <  1 ,  and  <V  1 , 


/3  4  a^/5  -  ... 


a  -  a?/ 


(43) 


Fence 


X 


(1  +  i  *rr)  [  (1  +  3*2  +  •..) 


5 


+  a 


3Wr 


£3 

'  1  +  Y 

(“- ) 

5k2 

-  ?^rtc2 

2  * 

—  Uj  Z  4  ... 

15 

'  V2  / 

2  ur  r  k. 

k1 

- 

+  i 


-  v2/* 

/</  4  “  k-  -f 


15 


5 

urc 


1  -  k 
k 


+  4  /4/T+  2k2 

h 


4 


(44) 


-12 

Now  for  metals  at  liquid  helium  temperatures  Z  >  10  “  second  and 

.  —^2 

0 ^  10  '  .  The  ultrasonic  frequencies  available  are  loss  than  10,# 
radians  per  second  and  hence  10  '  •  Thus  from  (44),  ^  1  and 

1,  all  the  terms  containing  factors  of  may  he  neglected. 

x  X. 


In  this  case,  k  V«lf  using  (28),  (41),  (40),  and  (44)  we  have. 


2k, 


m 


2  aZ 


l*t  (X^-  l)  *  m{*r 


/$  kXM Z 


kiM  1  « 


1  -  5k 


2 


2^rk 


Solving  for  od,  =  2k^  yields, 

x  _ 

^•1  4/  6.  / 

15M  V  3M  / 


-  / 


<.  -  *  4mvo^h-, 

j 


'45) 


°9 


O 

L 


since  k  ^  to-  and  V  ^  +  mv  as  is  seen  immediately, 

1  _  c  {-  S  _ O 

Je  5» 

Again  neglecting  the  terms  containing  ,  we  have 


—  m/S. 

K  *  '  < 


1  -  k 


2 


_2 

ki 


E^ 


rd 


X, 


'fi\ 


1  -  k 


2 


2 

-  v.ur 


Kk 


mvc  l1  -  4  \  +  4mvo  v-rkr+  2k2 


3M 


15M 


+ 


Therefore 


*/  1  +  m 


k? 


kJ 


l/2 

\ 


+  mv 


3M 


+  lu>'-zl  +  «r2* 2\  *  ••• 

15M  V  k. 


x  /  1  -  k 
k 


and 


V2 


3M 


to  obtain 


-  i 


Therefore  , 


A 


/i,  0  -«^/4k2) 

r 

1  4- 

2 

mv  + 

0 

4mv2  prW 

(l  +  m/l.T) 

15M  \ 

(1  -  <x'/4k2  ) 

'  1  + 

(v  f  + 

^  )  +  •  • 

/  (1  +  m/M) 

kl 

iinilarly,  for  k^/'>)l 

and  arbitrary  we 

*  ^  -  1  +  1 ,  + 

2  a  3? 

®  •  • 

P  1 
-  2k^k0  /  +  rt-. 

.  ,  2 
1  1  ti 

-  1\  +  ...  ■ 

Or  Z  6 

3  \  4 

/  J 

\  + 


-1 


(46) 


+  i 


# 


(47) 


-  dV2  / 1  - \  +  i>  k,  /  + 


3  4' 


k2 

'2 


3 


+ 


Hence  and  X .  v  k,  f  x  10  and  thus  we  may  again  neglect  the 

terms  containing  £k.  in  (40)  for  practical  frequencies  (ar(  10^  radians 

T 

per  second). 


and 


“ 


V 


Therefore,  (28),  (40),  and  (47)  give 

2 

.°  \ 


^  2-1 

//  bv  /  ^  +  mv  \  *  //  mv 

^  /  Z-' 


O' 

6M 


3M  j 


o. 


MV? 


/3M) 

1  - 

l 

°ct  +  •  •  • 

T 

t— 1 

1 

w* 

_ » 

‘l  +  mi?1-  l\ 

Ml  4  /. 

1 

w 

2k^ 

1 

£  1 
l—*  f\> 

'48) 


(49) 


Finally,  from  the  above  we  see  that  we  may  set 
2  2 


m 


m/z 


V 


-  1 


-  1 


\  g. 


1  -  (k^  )  1  tan_1(k1i’) 


(50) 


over  the  whole  frequencs'-  range  and  where 

g  «* 

We  note  that  (50)  is  precisely  the  result  obtained  by  Fippard,  ^ 1955 ) 
[6]  ,  and  Steinberg,  (1958&)  [  7  ]• 


The  development  is  just  the  same  as  that  used  in  the  last 


section  for  incident  plane  dilational  waves.  e  consider  a  plane  shear 
wave  travelling  in  the  x-direction  and  v/ith  a  displacement,  s,  end 


particle  velocity,  u,  associated  with  the  wave  in  the  y-direction.  In 


this  case  there  are  no  density  changes  and  hence  no  electric  fields 

resulting  from  a  space  charge.  However,  the  total  current  density, 

— >  — * 

j  <=  J  -  H  eu,  will,  in  general,  be  different  from  zero.  Tn  this 
e  o 

case,  magnetic  fields  will  be  generated  and  from  these  electric  fields 
by  induction,  mhe  effect  of  the  magnetic  fields  can  be  neglected  both 
in  the  Boltzmann  transport  equation  and  in  the  equation  determining 

the  ion  motion,  as  we  have  already  demonstrated,  Hence  we  have  only 

^ 

to  determine  £,  J  and  £  will  both  be  parallel  to  u  or  to  the 
y-axis .  Further,  div.£  *  0  and  the  equation  of  continuity.  (Pla), 
requires  that  £  »  0,  (no  space  charge). 

From  (19)  we  have 


(51) 


(1  +  iurt)  (1  -  ia'cose) 


where  vy  *  voinecor.^5  and  everything  else  has  been  defined  in  section 
2.2.  From  (?0), 


J  »  N  e  ez£"  +  u  1 
e  o 


(52) 


m 


where 


b  3  3  1  -  g  -  £ 

2(1  +  i-T)  [  a: 


# 


(53) 


2 


a 


TTence  , 


31 


G 


£  =  mu  k2  (l  -  b)  £  k2b  +  k2  -  , 


where  here,  k‘ 


4  Tfi 


and 


2 


N  e2r 

o 

m 


(54) 

(55) 


2  2  P 

ince  c»V  ,  we  may  neglects,  c  compared  to  k  throughout  the  rest  of 
this  section. 


Continuing, 


N  mu 
o 


(i 


-  b)k2  [k‘b  +  k2]  * 


(56) 


hun  the  equation  of  motion  for  the  ions,  (24)  become: 

#  . 

z38 


=  -  e  -  K  mu  (l  -  b)  (kc  4-  k‘  ) 

Pf  l.  °  o  v  y  v  o  *• 


(k2b  +  k2  ) 
v  o  ' 


(57) 


In  (57)  <5  is  the  square  of  the  velocity  of  the  s}iear  waves  in  the 
absence  of  the  electron-ion  interaction  and  k^  is  the  reciprocal  of 
the  classical  skin  0  epth  of  electromagnetic  waves  of  angular  frequency  ur. 

The  attenuation,  ,  and  the  velocity,  Vs ,  of  the  shear  waves 

depend  upon  the  relative  magnitude  of  k  and  kQ.  ‘  re  consider  rirst  the 

!2  2  I  2  2  | 

k  /k  I  and  jk  /k^b |  are  both  much  less  than  unity.  Now 


(57)  becomes 


y^ 


2  X 


i*^  mo 
_ cq  __ 

L- 


(58 


Setting  b  »  b^  -  ibf,then 


2k, 


i4/~n 


- 1 


,2  ,2 

b  +  b. 
0  1 


(59) 


3? 
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l+jr 


kg) 


'!3 


m  r 


,2 

b  + 
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1  J 


'60) 


There  are  now  two  cases  of  interest  in  this  region,  ^hese  are.  first, 

k^<<!  and,  second,  k/~l. 


For  ki<<l, 

b  =  1 


(l  +  ±i<s~z) 


Thu: 
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2  k, 


urtrik^f 


1  -  a 
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2  2 

=  rav^~  jr 
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i*r< 


/ 
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35 
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m  r  l 


2  2  2 
*rr  -  k,  ^  /^rT+  k 


P  +  2kp 
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(62) 


_  +  2«rV?k, 


-  8k4i4 
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■+* 


Therefore, 
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7*2 


Ard  -  kg/k^) 

(1  +  m/M) 


.  2  2  2 

1  +  mv  w  z 
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5M 
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+  <X  4- 


(1  +  m/M) 


2  2  2  2 

1  4-  rav  wr 4  oc '  4-  . . . 


>K 


r  63  > 


For  k^'-l,  then 


.1 


2(1  +  1  kt  r ) 


j£_-  2  4-1  I+^  + 

2  ~~2 

2a  a  &r  3  2 


(64) 


34 


Hence , 


oc 


w  in 


K^kj 


4k^/  +  ... 


;// 


4rnv  us- 
o 

JnUVJ 


7  65) 


and 


2 

ur  = 


(1  +  l6m/3/72M) 


1  +  m  f’4vQkri/ 


*u  t  > 

M  1 


3^VC 


2-16 


3/7 


+  mvo  Y2Q  {u?c7  -  1)  /  1  +  2  -  16 


1TV5  k^P 


3tt  jr ) 


y... 


(66) 


As  k-/  increases  further,  it  becomes  invalid  to  neglect  k 
compared  to  Jbk^  .  The  critical  frequency ,  ,  at  which  k  =  jkb  jis 

given  by 

%'L 


Ur, 


2  3 

3/r  v-V 


So, 


•  67) 


To  obtain  this  value  we  have  used  (64)  and  have  set  b^  =  J>7ft 

4k  f 

9 

For  most  metals  jfc  ^ 10  cycles  per  second.  At  frequencies  close  to 

277 

the  attenuation  starts  to  diminish  with  increase  of  frequency  until 

2  2 

eventually,  when  k  >2k^  and  k/T»1,  it  becomes  independent  of  frequency, 
cx.j  .and  Vs  are  now  given  by 


*s  * 
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(6fi) 


MVj  t 


and 


V2 
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1  +  3m'vf  1 

~2  2~2 
Mv  uf  £ 
o 


(69) 


35 


It  will  be  noticed  that  at  low  frequencies,  k^<l.  k^~l ,  and 
^ ^ « 1 ,  the  attenuation  and  dispersion  of  shear  waves  are  due  to  the 
induced  electric  force,  Whereas,  at  very  hi^h  frequencies.  */r>> 1,  they 
are  due  solely  to  the  collision  term,  I  ,  that  is,  they  are  due  to  the 

transfer  of  momentum  to  the  ions  from  the  electrons  by  direct  collisions. 

* 

This  is  in  contrast  to  the  case  of  dilational  waves  for  which  the  at¬ 
tenuation  and  dispersion  arise  from  the  induced  electric  fields  at  ell 
frequencies. 


2.4  DIoCUC  ION. 

The  expressions  for  the  attenuation,  (45)  and  (48)  for  dila- 
tional  waves  end  (62),  (65)  and  (68)  for  shear  waves,  are  just  the 
same  as  those  obtained  by  Pippard,  (1955)  [6],  s*nd  Steinberg.  (1988a) 

[7].  In  Pippard' s  case,  where  he  considers  the  relaxation  time  for 
the  electron  velocity,  it  seems  reasonable  to  identify  such  a  relaxa¬ 
tion  time  with  Z, ,  that  appropriate  for  electrical  conductivity,  how¬ 
ever,  as  pointed  out  by  Ha son  and  Pommel,  (1956)  [4].  if  this  is  done 
and  if  the  theoretical  attenuation  curve  for  tin  is  fitted  at  low 
frequencies,  ls.4  <<1  and  ufZ<<  1,  then  the  theoretical  curve  is  small  by 
a  factor  of  about  2.5  to  7.5  for  dilational  waves  in  the  region  k/  > 1  and 
it  is  small  by  a  factor  1.5  to  1.7  for  shear  waves  in  the  region  k/,  >1, 
urt^ 10 


Steinberg  has  given  an  explanation  for  this  discrepancy  by 
int'-rpreting  the  relaxation  time  appropriate  for  the  distribution  func¬ 
tion  and  therefore  that  appearing  in  his  expressions  for  a  ttenuation  to 
be  Zy/2.  If  we  take  expressions  (45)*  (48),  (62),  (65),  (68)  and  set 
Z  =  ZT,/ 2  and  follow  Mason  and  ’’dmmel's  procedure  of  fitting  the  theore¬ 
tical  curves  to  the  experimental  curves,  then  the  attenuation  will  be 
just  twice  as  large  as  that  given  from  Pippard' s  formulation  in  the 
higher  fre  uency  range,  k-/>l  and  urC,  1 ,  and  four  times  as  large  for 
shear  waves  at  extremely  high  frequencies,  kd,  >  >  1 ,  ^rs  >>1 #  However, 
fteinberg's  argument  is  erroneous  as  can  be  seen  readily  by  verifying 
that,  once  one  has  characterised  the  collision  integral  by  a  single 
relaxation  time,  as  in  (17)  f,  v  and  all  relax  toward  their  equili¬ 
brium  values  at  the  same  rate0 


7  f 

J 


7 


We  shall  3how  in  the  following  sections  that  the  explanation 


of  the  experimental  results  lies  in  the  fact  that  the  collision  in¬ 
tegral  cannot  be  characterized  by  a  single  relaxation  time  for  the  whole 
frequency  range.  In  particular,  it  turns  out  for  the  low  frequency  re¬ 
gion,  k/  <<  1,  <<  1,  that  a  relaxation  time,  Zz  ,  appears  in  the  ex¬ 

pressions  for  the  attenuation  and.  we  may  write  7^  rat,  where  £  lies 
between  1/3  and  1,  The  precise  value  of  £  depends  upon  the  details  of 
the  scattering  mechanisms, 

3efore  we  teke  up  the  calculations  of  the  relaxation  times,  we 
mention  that  our  results  for  the  dependence  of  the  velocities  of  the 
sound  waves  on  the  frequency  arc  qualitatively  the  same  as  the  results 
obtained  by  Steinberg.  Since  the  velocity  changes  over  the  entire  fre¬ 
quency  range  are  of  the  order  of  the  mass  of  the  electron  over  the  mass 

-5 

of  the  ion,  that  is,  of  the  order  of  10  ,  we  shall  not  discuss  them 

further.  We  should  note,  however,  in  comparing  (46)  and  (49)  with  the 
corresponding  expressions  in  section  4  of  Steinberg's  paper,  [7],  one 


since  this  is  the  effective 


must 


elastic  constant  seen  by  the  dilational  sound  wave.  It  is  interesting 
to  note  that,  in  the  case  of  the  dilational  waves,  the  electron-ion 
interaction  contributes  a  term 


to  the  velocity  independently  of  the  r elaxation  processes  due  to  col¬ 
lisions.  This  term  has  its  origin  in  the  fact  that  when  a  dilrtional 


wave  passes  through  a  metal  both  the  ionic  and  electronic  densities 


vary  from  point  to  point.  Since  the  electrons  obey  the  ^auli  exclusion 
principle,  their  contribution  to  the  bulk  modulus  becomes  appreciable. 
In  the  case  of  shear  waves,  no  density  variations  occur  and  hence  no 
equivalent  term  is  obtained. 

"  If  we  do  not  neglect  the  term  appearing  in  the  deno¬ 

minator  of  (39)  but  do  neglect  the  terms  which  give  rise  to  the  atte¬ 
nuation,  we  obtain 

v  a  /?  +  mv  (1  +  *  v  k  /  3 ) 

y  rte  _ o  v  o  /  ^ ' 

3M 

This  expression  is  approximately  valid  at  all  frequencies  as  long  as 
the  metal  can  be  described  by  the  Sommer f eld  free  electron  model.  The 
above  expression  has  been  derived  by  3hatia,  (1955)  [lo] ,  in  a  much 
simpler  manner  in  another  connection. 
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.  CALCULATION  OF  TT-;  RELAXATION  TIMES  FOR  THE  ELKCTROn-FTTHON 
INT  .TRACTION. 

3.1  FORI-.IULATI  ON . 

We  again  formulate  the  problem  in  terms  of  the  W:  xrwell- Boltz¬ 
mann  transport  equation  find  employ  the  same  physical  concepts  as  used 
in  section  2,  Instead  of  making  the  assumption  that  the  collision  in¬ 
tegral  is  characterized  by  a  single  relaxation  time,  v.e  set  ur  the  col¬ 
lision  integral  so  that  it  takes  into  account  the  collisions  between 
the  electrons  and  the  lattice  vibrations.  To  do  this  we  start  with 
equation  9.33.6  in  Wilson's  "Theory  of  Metals",  (1954)  j  11 J „  and  pro- 
•ceed  from  there.  Using  the  assumption  that  the  equilibrium  distribution 
in  given  by  (13)  we  set  the  collision  integral  zero  for  this  function. 
The  experimental  data  also  suggests  that  there  is  very  little  energy 
dissipation  due  to  thermal  currents  and  hence  we  use  the  phonon  distri¬ 
bution  which  corresponds  to  thermal  equilibrium  in  the  undisturbed 

(Se ii) 

state  of  the  metal/.  Equation  9.33.6  from  Til  son  reads 


~}f 

0  C2  c- 

1 

[{ 

at. 

C-< 

(?G  +  D3  ®  f 

4-  v4/2 

c* 

^(f)(l-f(f^)]A(y^r  +h>J.)  +  |Njf(K-q)(l-f(K)) 

-(N-+l)f(KXl-f (X-q))f  j  , 

where  the  are  unit  vectors  defining  the  direction  of  vibration 
for  each  of  the  three  acoustical  modes  designated  by  j.  mhe  3(s-l) 
optical  modes  are  absent  in  a  metal  of  all  the  same  type  atom.  ~q 
is  the  phonon  wave  vector,  M  in  the  mass  of  the  ion.  h  is  Planck's 
constant,  V+j  is  the  frequency  of  this  lattice  vibration,  (2G  4-  l)^ 


(70) 
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is  the  number  of  atoms  in  a  cyclical  three  dimensional  crystal  forming 

2 

a  Ir avals  lattice,  and  C‘  is  a  constant  measuring  the  strength  of  the 
electron-ion  interaction  and  is  of  the  order  of  a  few  electron  volts. 

The  summation  goes  over  all  possible  values  q  to  take  into  .account 
both  the  emission  and  the  absorption  of  thermal  phonons  by  the  electrons. 
The  represent  the  probability  that  there  is  a  phonon  designated  by  q 
at  a  given  temperature,  T.  Tor  thermal  equilibrium 

-  [exp(  h vf/kT  )  -  l]"1,  (71) 

where  k  is  the  Boltzmann  constant  and 
o 

_fl(x)=  sin(  xt /#  )  ("J2) 

x/Vr 

and  hence  acts  as  a  delta  function  conserving  energy,  f (ic) 
is  the  electron  distribution  function  and  represents  the  probability 
that  the  electron  state  designated  by  K  is  occupied.  ',Te  note  here  that 
there  is  a  two-fold  spin  degeneraoy  which  must  be  taken  into  account  but 
vill  not  appear  explicitly  in  the  soltzmann  equation. 

If  we  now  set 


f 


o  ? 


>>  E 


and  set  the  collision  integral  equal  to  zero  for  f,  then  (70)  reduces  to 
9.33.7  In  \'il son,  that  is. 


If 
bt 

Co U 
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xA(E.-E,  .+hvO  +  f  (K+q)  \  1  -  f  (K)[ *-hvf) 
'a  A'*  s'  o  '  C  o  J  K  Kr'i 


f2(5c+q)uifo(i';A.?)  -  f2(ic)«f 


O 


^  C  /<  r  1 


d£ 


dq 


(74) 


where  we  now  use  f  rather  than  f,  since  we  are  interested  only  in  a 
first  order  theory  in  u ,  and  A  is  the  volume  of  a  unit  cell. 

In  working  further  with  this  integral  we  shall  assume  that 
Er*  is  a  quadratic  function  of  j  PC  j  and  has  spherical  symmetry,  and 
that  is  o  function  of  j  q  |  alone  with  the  values  of  |  <f  |  lying  inside 


a  sphere  of  radius  q^ .  For  an  elastic  continuum, 


o 


(6H2)/j  , 


q 


a 


% V--  <  hvo 


y.o& , 


(75) 


where  js  the  maximum  frequency  supported  by  the  lattice, @  is  the 
Debye  temperature  for  the  solid,  and  a  is  the  linear  dimension  of  a 
unit  cell  in  the  solid. 

We  now  write  q  in  terms  of  where  o  is  the  length  of 

q,  c  is  the  angle  bctv/een  q  and  V,  and  a/-  represents  an  azimuthal  angle 

— ¥ 

about  K  taken  as  the  polar  axis,  "ith  this  in  mind  we  can  writo 
Kx  <=  PCcos  e>  j 


X 


y 


K 


and 


q , 


t>  JCsin<9cos/5  ) 

»  Kainesi n/  ? 

q(cosaco3^  +  sin^sin^cos^)  ? 

*  q(sin#cos/costf'  -  cosecos/sin/cos^  +  sin^sin^sin-^)  > 

*  qr sinesin/cos y  -  cosesin^sinycos^  -  coe/sin2(sin<*')  « 

:.'inco  the  A’ s  act  as  delta  functions  inmergy.  then  we  can 
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carry  out  the  integration  over  4  and  we  writs 
cos^  =  (  h>t  -  E-)/*v q  , 

and 


con  *  (  hv^*  +  E^-)/£vq  >  ^76) 

where  the  plus  sign  on  £  indicates  the  absorption  of  a  phonon  by  the 


O  o 


election  and  the  minus  sign  the  emission  of  a  phonon,  ft-  =.  -bo' 

*  “2m 

the  energy  of  an  electron  with  the  wave  vector  q. 

As  demonstrated  by  Wilson,  we  must  set 

E-?  h  >£)  dco&ft  «  n(qddcosy  +  Ea  +  h2.)dcosft 
*  1  (  y  k  ^  d ) 


h 

2Vq(dE)" 


3iny  dy  -  h  S(y)  dy  9 

y 


^  c:IT 


where 


and 


K 


=  2!7t  (  E  +  hvi  +  oddcosft  ) 

T*  ^  t7 


2  2 

=  *l  K  - 
2m 


In  this  event, 


af  ■ 

»  C?A  r o/dqda,- 

'(Hs+  1)  (e?  +1)  £ 

.at 

v  ■  /  . 

87*  ?5hv  ° 

(o’-  +  X)  " 

f2(?) 


-  q  (K+5 


uif  +  K 
o  5 


d  E 


is 


(77) 


(e*  +1)  {  f2(K)  -  f2(K-q_)i  uJf,' 

(o’'*  +1)  ^  35  . 

where  the  plus  sign  indicates  that  cos ft  must  bo  used  and  the  minus 

T 

sign  that  cooft.  must  be  used.  "-Iso  we  have  set 

and  z  »  hv4  <• 


■''78'' 


(  S  -  .30 ) 


(79) 


k  T 
o 


k  T 
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Restricting  the  values  of  q  as  in  (75)  we  can  reduce  r78)  to 


df 


3t 

co/t 


(«?  +1)  !  f,(K) 


(80) 


where 


1 

A 


o  p 

3h^C 


4) 

\4//  /  4Hak  $ 


and  the  last  term  indicates  the  average  value  over  ur.  In  setting  the 
limits  of  integration  wo  have  used  the  fact  that  for  metals  ?T>  Nc/4 
and  for  this  case  the  restriction  on  v  is  k  However,  for  semi- 

metals  and  semi-conductors  IT  <  17  /4  and  the  more  restrictive  condition 
hv^hv  K/tt  must  be  used.  This  results  in  the  limit  (S)  Ap  being  replaced 

o  '  ‘  j 

by  l2]$i/l(£)  nd  this  generalization  can  be  readily  carried  out  in  the 

\  DJ  T 

subsequent  work,  although,  it  h?:s  been  omitted  here.  are  now  in  a 
position  to  start  our  solution. 

j^As  v/e  have  mentioned,  Holstein,  (1959)  \  12 J  ,  has  very  recently 
proven  that  f,  (l4)*.is  actually  the  distribution  function,  corres¬ 
ponding  to  local  equilibrium,  towards  which  the  electrons  relax  in  the 
presence  of  the  ultrasonic  wave,  through  the  electron-phonon  interaction. 
Before  proceeding  further  we  shov;  that  Holstein’s  collision  integral 
reduces  to  (80)  above  on  which  our  considerations  are  based. 

He  finds  that  one  must  szt  the  lattice  displacement  vector 

R->  equal  to  the  sum  of  the  displacement  due  to  the  thermal  motion, 
a 

R~^  ,  and  the  displacement  due  to  the  ultrasonic  wave,  s,  that  is. 


V  rs. 


R  -  4-  s .  He  then  considers  the  change  in  potential  energy 

0  ■' 
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isplaoem  it  >f  th  ioni  y  in  ti  n  ■.  ?-r 

in  t  ■  r  3eri<  »  in  H*  ^  ana  s  ma  retainm  ■  1 1  lm^;r  in 

9  /r«  J 

'  io’  is  of  t;  order  of  the  theory  c  mnidered  ’’.ore.  The  o’  r  ••  .  in 
tho  potential  energy  cont  i  the  u  u.il  thermal  terra  linear  in 
an  addition.  1  t  rr  t  correspond ing  to  V  o  eirroilianeon:;  Imorption  of 


a  thee  al  a  aoui  ion* 


i  '  jur.t  Vi:  addition  al 


vm  hich  c)  n  t  cmilibrium  listribution  fr-  .  f  to  f.  i:  'inde 

y 

t  at  one  m  ,  us ?  just  (70)  to  determine  the  colli  ion  in  ;gral  >xc  *pt 
that  now  the  anergy  conservation  terras  mu  1  b<  raodifj  to  re 


a{'  K^-  :s  -  hv/4  +  «*•(**-  Vj^/  )]  Szti<’-$-K(ZZ') 


and 


id  [  '  -  -  -  ^  +  '  1 

<■  /<~q  k  1 


+  rrai  •  (  v  —  —  v  ^ 

mu.  \  ,  t<  f< 


y  ,  oit>) 


./here 


•v  f T'  ?') 


(v-  -  V-,  <?.  r0) 


•  nic) 


and 


t  q  ■ ' 

"e  .  (r  }  *  i^jun*.  e  ‘  +  complex  conjugate 

<-  v  O  '  *  4 


^  31  -  •  'I 


o*  (r_)  i  tl  5  tr  in  d;  lie  ■  ■  ci  ted 


L  ^  n 

and  is  pronorti'  nr.  1  tr>  the  lattice  dia;  V  c  >en t . 


ubotituting  ^14} » 

n  |  v  -  v\ 


[1(? 


5, 


1/vi  +1 


(81a)  to  (Bid),  ond  (71)  into  (70),  we  see  tint  f  i-.  tie  cnilis.riuu 
distribution  function,  und  r  tbo  aboV‘>  mosi: ication  o  •»  .  on.  ?  .'r 
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the  .first  terra  in  (?0)  subject  to  t 


»e  ncv;  energy  ro:t?5.ntionr 


ks  •  [i  -  f(")j  -  "af (")  [  i  -  r(v+r)j 

s,+  hv,»J'V) 


p^t  'Y 
x  z1  W+f  / 


/ 

1  1 

2 

7 

1 

0 

1  ■*/ 

.  2 

fe1 

v  - 

u  |  - 

-  f  /k 
o  J  »  /  /  o 


*  a  f  (*  )f  (‘  4f)  ;xpf/ m  | v  -  u  |  2  - 


i,  +  hh 


r/  ut7 

■*  |  2  1  -* 

— ^  , 

r/m|v 

-  u  —  m  v 

-  u 

/  1  — 

1 

t\2 

2 

*=  o, 


luce 


,<  +  h  ^ 


nru  *  v-. 
'  A 


-  v 


/<'  > 


inh  so 


m  I  v 


u 


-  in  v  -  u 


-  hy> 


-  h  vt 


)  /  V  } 

)  ^ 


k  v> 

q 


+  DIU .  r  V, 


l< 


-  V' 


=  ..o  , 


nct't  th ;■  first  tern  in  (70)  is  zero.  imilsrly,  t ’•<?  second  :  rm  :: 
loo  s.:ro  for  f  =  f. 

(73)  in  (70)  an''  ”ol  toir.’  result  ,  th. 


t'inally,  if  /e  use 


:olli.  ion  integral  reduces  to  (74)  with  f  repl  ced  by 


■  JV 


nf 


r  •  lr  ced  by  (81a)  rm  (01b).  Yirther,  remembering  tl  st  w  :re  r  t.  ining 
only  terms  linear  in  u  and  noting  that  all  the  correction  t  \  re  ro- 
portionol  to  u,  f'  have  to  first  orde  *  in  u  just  t’  o  o.-cyr --  ion  /'7T,>  r 
it  stands,  so  that  the  collision  integral  is  given  by  'SO)"]. 
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Proceeding  now,  it  is  more  convenient  to  take 

f  =  f  -  f.udf 

o  1 _ o  , 

3E 


on  the  left-hand  side  of  the  Boltzmann,  equation. 
(73)  with  (82)  v,-o  h?<ve 


Ci  mrnrin* 


fiu 


H>  — 

rnv  #u 


^l+  ''2U 


ubstituting  (82 )  and  (03)  into  (80)  results  in 

2 
7 
It 

Ti-7^ 


3f 

3t 


V 


dz 


(C*  4-1) 


-  / 


ra  ) 


m 


' p 


mv.u  IscosY  \+  f ,  (")u 


x-  Y- 


=  -mv.u  ^  f 


v;  here 


+  u  3f  I, 

o  1 


Z'/v)  3 


say, 


1-c"2  |  (e^*  +1) 


'<v  > 


(£_+!)  ffjO )  -</•  " > 


^trid 


(g?  +1)  f  1 


-2‘  ^f%i;i2: 


^  •>  J 


__  V 

K  » 
O 


'8’"b> 


/here 


4ma2 


Wo  note  that  (85*0  in  identical  to  the  ex'  re  or-:  ion  found  for  f  r  - 
tion  tinii  fo  el  otrical  conductivity  on  rewriting  9«  ;  '  e  f  Wileon,  [l2J 
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'he  complete  loltzmann  equation  can  be  Tit  -n,  us  in  '32) 
/\nd  (84)  in  16),  in  the  fox-m 


-i uj  +  iVv>,)f/u  3f ,  +  et.v  *  -  inv .v  3f 

3  3  3 E  3  1 

'e  2iov/  expand  f  in  terms  of  spherical  harmonics, 

tinu  jus  function  this  is  nV  ays  pos  Lbl<  and  a 

such  9  well  behaved  solution,  ’once, 


(8f) 


i7'  >r  n  hounded  eon- 
•unv  thr  t  '06)  has 


•u 


-  n 

-f.  m 


(*,/) 


'37) 


^he  next  step  is  to  substitute  (87)  into  (°6),  multiply  throu  ■  by 
Yp  (<Sq^)  •  ineded^  ,  and  then  integrate  over  all  male  , 

is  gives  us  the  fellow  in,;,  integrals  to  ev'  lu  L  ^ : 


/ 


-g  rn 


Y/  ( e *0)  (v)jJ,(&$/)<p±ned&d/ 

JT&n  (p+n )/  v?0*(v)  ; 

(2p+l)  (p-n; / 


2 


where  ^  is  one  for  n^O  and  two  for  n  =  0 . 


jV-x (@,/4)sined©d^ 

^  n 

J  v  (<?,^)sin<?ded^ 


2  TT  4* 

T^pTi) 

2  Aj 

T2ST) 


tn+n)-  v  d-  ho  • 

(P-n).' 


-v  V  S 


JE2L: 

(p-rnZ 


r' 


/?  / 


f  yxYp(e*rt  5Z  v^cJT  (v)Y^  (e^)*in&ted/ 


KW)  i_  v'"c"(v) 

(/ +m)Y^_  , 

(®.y) 

+  (^-r-1 

/ 

j 

X 

r;ine?r  6»  <jt>  , 

f  o 


8a) 


'38b) 

'88c) 
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-  2  &  a,-,  ( p-t-n )  l 

(2p+l)  (p -n)' 


i  \  ?  /~i  ^ 

JlEi*LLv  V  p_, 

( 2t)-l) 

L  \  / 


fr2  n 


(v)  +  (p+n+1 )y'  ~ c 

(2r>+7,) 


P+t 


(v) 


J  i 


'op.;  \ 


V,rts?r,  in  fcl  oolli  ion  integral  m  no  c  n-y  ut  ’  2 
intogn  tion  ov  r  :  nd  tlir  ,  ;vos 


£ —  ^  •  \L'P  9  ' '  u. r 


7”  ▼ '  '<£(▼ ' )  p  (cor  A  )T%,/) 


'Y’y  /X 

""ho  steps  in  (88e)  follow  ince  v/,  =  vfl  -r  3: 
and 

*  y(cor,n)'7(e^) 


'38o) 


fc. 


is  inch  p';nd  nt  of 


'p.flf) 


Cn  these  expressions  n  is  now  the  angle  betwe  in  k  =  ]:  +  3? 


In  this  case 


con  fl  * 

f  1  +  k  ff’z  -  D/'  \  2 

r  1  +  k  "yz  1 

0  1 

0 

2B  E^/ 

E 

- 

1  '2 


v 


v 


1  I  V  r"W 

X  T  i'.  < 


t2  ?i 

z 


2-; 


'83t-N 


Substituting  (88a'  to  (08g)  into  (80)  on  i'i  in;  out  th 

# 

faotor  6m  (jo _ i  n)  1  the  oltzmann  equation  becc 

W+l)  tr  -  n;/ 


-  iurv^c”  (v)  +  ik 


r  f 


s  P  n 


LfcIUv  V'  (v) 


(p+n+l)v',"20  ”  O'! 


(2]~1) 


(2p+3) 


/>-/ 
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+  ev 

V) 

i 

+ 

to 

-  O'  V/  ^  / 

+  AH  1 

pi  n-, 

=  -  nv 

u 

'  J 

A  . 

£  S  + 

°ft  no 


/T7) 


s/>,i,  +  VV' 


.,i3 


<^/r 


2 , 

z  da 


(e*+l) 


A^i  v  \®J*y£yT  |  1-e  ^  |  (e^^+l) 


A 'n  0  /  N 
n  IV 


P 


,P  « 


v'  ^(v')^(cobA) 


t  o 


9) 


can  now  notice  several  features  of  the  ^oltsm-nn  equation 
in  this  form.  ]ach  p  and  n  value  gives  a  linear  inte  :r  1  equation 
coupling  three  consecutive  C^’s  with  the  same  n.  Indeed  e  ’  ave  ;?n 
infinite  set  of  coupled  linear  integral  equations  for  each  n.  Frrtf  r 
if  we  impose  the  conditions  of  no  y~  and  ^-components  o'.‘  cusran^  "or  the 
incident  plane  dilational  wave  moving  in  the  x-direction,  that  ir  C;  - 


- 1 


<  =  0,  then  we  ’nave  all 

Cj1  =  0 

r 

hence  in  t>  is  case  we  need  only  retain  the  set  of  equation;'  "or  n  =  0. 

imilarly,  for  a  plane  shear  wave  moving  in  the  x-direction  n  pol  risen 
in  the  y— direction  we  need  only  consider  the  set  of  equation''  for  n  =  1. 
bus  we  have:  (i)  for  the  incident  plane  dilation  !  waver: ; 


frti 

i. 


-  wCp'v)  +  ik 


n  2, 


P  <v_,(v)  +  (p+l;V  (v) 


=  -  m 


2t>~1 


ti\3 


(2p+3) 


tp 

+  c£ 


$ 


ft 


u 


®/j 


2Zd3 


(e^+l) 


r.  fy)  - 
p  ’ 


/\fiZv  \ ®)  _®/j-  |  l-o  |  (e^^+1; 


vY  Vp  (cosfL)Cp(v/) 


90a} 


and 


where  p  =  0,  1,  2,  3>  •••  » 


(ii)  for  the  incident  plane  sheer  waves; 
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-  i««rC/>(v)  +  ik  T  (  r'~l ^  C  (v)  +  (p+2)v*CL  ^v) 

'  TS&t 


eg 

u 


pt 


®/r 


*  -  n  §  +  _JL _ / ; 

r.  "  A ^w/Jq. 


zZdz 


1-e 

»» 


-z 


(chi'!  f 

(e?«+l)  L 


<vM 


vM  ^p(coeil)C/?(v/) 


'90b) 


where  p  =  1,2,3,...  ,  and  the  index  n  =  1  has  been  suppressed  in  (90b). 

To  deal  with  the  integrals  appearing  on  the  right  bond  sides 
of  (90a)  and  (90b)  we  assume  that  each  Ch>  (v)  is  a  very  slovly  varying 
function  of  energy  so  that,  on  integrating  these  expressions  multiplied 
by  Zf  v'dv,  over  all  energies,  we  can  set  CP(v')  2=r  0  1  (v)  in  those  in¬ 


tegrals.  We  are  then  left  with  a  temper,  tu re  series  tl  rough  the 

.  p  ,  ^ 

v '  x  Pp  [cos  12 ;  and  since  e  are  considering  only  v  -ry  low  tem- 


l 


p< ratures  we  need  only  retain  the  lowest  order  temperature  terns,  ^be 

next  step  is  to  multiply  (90a)  end  (90b)  by  v  hf  dv  and  integrrte  over 

o 

3E 

all  energies.  r”his  integration  is  considerably  sin  lifiod  by  the  fact 
that  df 


3 


is  almost  an  energy  delta  function. 


‘'e  have  now  to  evaluate  integrals  of  the  tyre: 


J  v'^cpvjjfgdv 


CC 


v 71  r  /  C  (v)5(E~E  )  dE 

P  v  '  o' 


d?  m 


v  ^  * '  C  ^  v  '>  , 
o  />'  o'  > 

&/T 
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mv  df  dv 

_  _ o 

-/  <>E 


2, 

v  dv 


m 


d* 


z  dz 


(•*+!) 


A-fiVTSv  dE  -*>r  |l-o'z|  ((/'z+l) 


'91a) 
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-  Vz 


Vz 


2r,/?\2z2 

01 


\S  rr 

A\.  CJ 

O 


D 


A  -f^E 


oc 


’si® 


(Olb  > 


n* 


T  7. 

\  af  V*  dv 
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(?■ 


W/r 


2, 

z  dz 


(e?+l) 


o  /\fi 1  v  '  A'  a  E  -  #V 


1-e"2  |  l  e' 


f-q+z 


•i) 


C  M 
P  ' 


/,  P 


x  (co  n^O') 
e?  d;p 


C/v) 


-  /v 


AP 


A  |l-c~Z  I  ,(e*+l)  (efrz+l) 

P/?  (cosxi)^(y/) 

'  \P 


a(v  )  5 1 

p  o'  * 


v  V  r«  / cor  fl v 
o  1  P 


1  -  ^(cosnJJl  +  kQTa 


O 


(91o ) 
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'"he  subscript  aero  denotes  that  the  velocity  and  ener  y  are  to  bo  eva¬ 
luated  at  the  Fermi  surface.  With  the  above  approximations  we  have 
reduced  the  initial  set  of  coupled  linear  integral  equations  to  simply 
a  set  of  linear  algebraic  equations,  in  the  C^fvJ,  which  we  may  now 
solve.  (91h)  follows  from  the  fact  that  only  even  power."  of  s  in  the 
curly  brackets  of  the  integrand  give  non-zero  contributions  to  the  in¬ 
tegral  and  from  the  definition  of  Jnf 
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®/r 

zn  dz 


We  notice  here  that  (96b)  is  simply 


where  Zt  is  the  relaxation  time,  due  to  the  electron-phonon  int  r- 
action,  which  appears  in  the  electrical  conductivity  at  lo  temp  r- — 

tures . 

"Ur t her ,  the  expression  in  the  square  bracket  of  foie)  i.> 
a  series  in  powers  of  z  and  T.  First,  as  we  hnvc  noted,  wo  ne*-d  re¬ 
tain  only  the  even  powers  of  z  in  this  polynomial  ns  t’  e  od  ’  powers 
integrate  out  to  zero,  secondly,  since  we  are  cons if  rir  y  very  low 
temperatures,  v;e  need  ret  - in  only  the  lowest  order  term  in  th  ■  tem¬ 
perature.  "it},  thi3  in  mind,  solutions  for  the  two  ones  or  ler 
frequencies,  l:P,4<  1,  and  of  ‘requencies,  k^>l,  obt  inod. 

As  a  preliminary  to  sections  and  5.4  0  mrke  tl  e  pol- 


lowing  id en  t i t i es : 
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?o  obtain  a  complete  picture  of  the  collision  integral 
an  calculated  here,  we  rewrite  (15)  in  the  form, 
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a  0  ,  for  p  =  2,3,... 

(82)  may  also  be  rewritten  in  the  form 
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^he  collision  integral,  (80) ,  is  equivalent  to 
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'93) 


only  if  7/e  male  the  identity 
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'  (v  )u 

o  o' 


(94) 


This  identity  is  the  equivalent  of  (2?),  that  is  the  equation  of  con¬ 
tinuity . 


We  note  that,  in  reality,  it  was  necessary  to  carry  out  the 
integration  over  energy  before  a  reduction  of  the  collision  integr/  1 
into  the  form  of  (9?)  occurred,  so  that  (95)  is  oertainly  not  rigou- 
rcus*  however,  (93)  is  equivalent  to  the  preceding  procedure  where 

et  dff  *  —  o(E  -  E  )  and  ta)  the  )p (v) ’ s  as  slowly  varying  funo- 

Tr, 

tion  of  energy.  We  also  note  that  0o(v)  need  not  he  a  constant  hut 
need  only  satisfy  (94)*  Hence  in  thin  case  we  may  consider  ore)  term 
in  the  expansion  to  relax  towards  its  equilibrium  value  ith  r  charac¬ 
teristic  relaxation  tine,  rp9  and  that  there  is  no  sir  le  re]  xation 
time  suitable  to  characterize  the  collision  integral. 

We  shall  now  consider  Tt  and  Zz  in  some  detail.  c  Icu- 

lation  of  Zt  has  been  accomplished  in  (91b)  and  is 
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In  the  same  manner, 
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fir  1  . 

rn he  second  step  follov/s  since  all  odd  powers  z  in  tie 
integrate  to  zero.  r'he  final  reduction  in  (96)  arises 
that  we  are  considering  very  low  temperatures,  o  t) . 

<  0.05  or  T  <  20°K  . 

<8> 


'QO 


square  bra  ckots 
from  the  f  ct 
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In  metals,  for  the  number  of  free  electrons  approximately  equal  to  t’ 

number  of  atoms  and  where  each  unit  cell  contains  only  one  or  too  toms 

§ 

then  1  os  E  .  Y>r  these  conditions, 
o 

(ko<0  )Z  „  1CT® 


since  (&  ~  300°K.  Also, 


>  57 


e  remark  that  expression  (9*5)  fo*4  ^  in  the  .ame  that 
obtained  by  Steinberg,  (1950b)  [9],  in  Ids  calculation  of  the  visco¬ 
sity  of  an  electron  gas  resulting  from  the  electron-phonon  interaction. 

Lext ,  we  shall  consider  the  relaxation  times  in  eni  rel. 

I,  is  defined  in  (91c)  and  is 
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where  P.  is  the  pth  Legendre  cl.nomial  and 
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where  w  may  no  v  take 
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=  1  -  $  ,  nay . 

We  may  now  write  !\,  (cosfl„)  as  a  po  er  serio?  in  S  r. m'  obt  in 
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'e  now  r  tain  only  the  linear  term  In  &  since  the  other  terms  cent?  in 
higher  powers  of  T  and  for  very  low  temperatures  these  ore  expected 


to  be  docrea singly  small, 


hence 
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^or  p  »  1  and  2,  (97)  reduces,  of  course,  to  (95)  and  (96). 

It  in  not  immediately  clear  that  .ve  are  justified  in  neglect 

in/;  the  higher  order  temperature  ter me  in  each  r .  We  note  how  ver 

^  » 

that  these  t  jrms  will  bo  of  the  order 


n  -y 


n  =  7,9,11, . . . 


jw©7?) 

and  we  have  evaluated  this  quantity  in  table  I  for  m  $ 


0.05 


'he 


last  row  indicates  the  order  of  magnitude  of  the  neglected  terms  rnd 
we  see  that  for  all  practical  purposes  we  need  only  r  twin  the  hirst 
term,  t  at  is,  for  T.£  20°K*  The  values  for  J„  {(£) /t )  ■  r  obi  ined 

from  Wilson,  [12J  (p337)» 

We  note  that  the  relaxation  times  for  higher  order  harmonic: 
in  the  expansion  bee  me  rapidly  smaller  than  that  for  the  first  term. 


T  ence,  at  low  ultrasonic  frequencies,  we  expect  the  higher  order  h:  r- 


mondes  to  be  able  to  return  to  t.'  oir  c  ui librium  value  and  thus  only 
the  first  few  terms  will  be  appreciably  different  from  soro.  Tn  the 
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next  tv/o  section.",  we  consider  this  case*  At  Mg’  or  Areouenci es,  the 
higher  order  terms  do  not  have  time  to  reach  their  equilibrium  values 
and  hence  they  must  be  taken  into  account. 
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3.3  lo  ’  ^3rTKKCY  I'  c JT)  :  i !. .  •:  r.r-avi  wwss. 

r.inee  the  shear  waves  present  somewhat  simpler  expressions 
than  do  the  dilational  waves,  we  shall  treat  this  ca  e  first.  for 
low  frequencies,  ki* < <  1 ,  we  make  the  approximation  that  Mr  a  ries  ox_ 
pan a ion,  (07),  is  rapidly  convergent  and  retain  only  th  first  two 
terms.  \s  the  solution  unfolds,  we  see  that  this  is  r  relatively  good, 
approximation  in  this  frequency  region.  Rewriting  r90b)  with  th  old 
of  (91a),  (91"b),  (91c),  and  (92)  we  have 
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where  ./  =  v  r  and 
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^  0  .  ^ 
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The  solution  now  proceeds  by  neglecting  !?/  compare’  to  J^/e 
in  (100).  The  resulting  equations  can  now  be  solved  for  J0/e  and 
and  yield 
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wl  ere 
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'e  nee  that  M^/j  ^  ik/  <<  1.  TVds  process  continuer  ro  th  t 

G  * 

G 

■•^KY,  1  nd  hence  we  map  neglect  r';  compared  to  J  /e  in  '100). 

The  solution  riov  follows  exactly  ti  e  ame  steps  :  n  those 
in  section^ 2  .  5  0  Thun  we  br-ve 
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(analogous  to  equation  (54)), 
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(analogous  to  equation  (56)),  and 
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(pnalogous  to  equation  (57)) • 

Since  we  are  considering  lo  frequencies,  t'  n  0  r  in  t1  o 
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the  ions  can  be  written, 
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'"o  obtain  higher  terms  we  must,  of  cour  e,  carry  ‘Vrther  t  rr> 
series  exp.'  nnion  for  f,  .  'g;  in,  as  previ  -u  l.:  ,  Cot'  V  #  <  <b  1 


in  t’  o 


1-1 


/ 


>5  ^ 


hkf  -  4)  (1  -  2^r2).(iv  ...  j 


63 


- 1  f?k,k* +  1  -  a*+  •••/ 


Therefore 


^  = 
3 


2h  ■  ^  yy* 

Mr  5k, ^  L 


fk 


kf)  n 

1  '  ' 


) 


2tsr  y  }  + 

z  /  z  *  *  * 


•ence 


ciVo  ujTz 

r*/* k 


k. 


k, )  , 


«  O  O 


cxh  -  nnrk.^j^ 


51' 


jn-tr  ■£  I  ^  —— 

UiiiLi ? 

5'  T-* 

5 


■10?) 


The  dispersion  relations  are 
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We  remark  that  (10?)  and.  (103)  are  precisely  th  •  or-  *  r, 

(6?)  and  (63),  section  ?o3»  with  Z*  merely  replaced  by  Z^0  '  t  i  , 

these  expres.-.-.ione  do  not  contain  Tt  ,  the  relax/- 1 ion  time  i- to 

for  electrical  conductivity.  Of  course,  for  lo  •  temperatures 


the  electron— phonon  interaction  alone  >ve  have  seen  t’ •  t  lz  i  very 

i?rly  JL  %  ,  in  gen  ral,  however,^  -  ill  not  he  uc!  a  ‘  function 
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5.4  J~w  FHBQVirCV  I'  :j’>  '  r’  '"LAN  fib A 7 I Of Ah  WAV  'S. 

Tere  we  deal  wit!  tlie  more  diffici.lt  case  of  0 ilatiom  1 
i  jg.  rhe  difficulty  arises  due  tc  tl  e  fact  that  V  bulk  r odulu* 
of  the  electron  gas  nust  be  token  into  account,  lie  ltir(  '90s  ^  wit! 
tl#.  aid  of  (91a),  (91b),  rnd  (91c)  vie  Lave 


(1  +  i^5)C/O(v0)  -  ikr 

=  n>(l  +  !Llr.)i 

\  mi*.  /  i  ) 

'or  p  =  0  we  must  use 


TT 


V) 

— 


IT 


p 


f\  )  +  f"  +  ljV^-C 

,v  o'  ^r~^T  . 


>> 


2,  +  ’T  ^ 


_  I  0  7 

~  -*->  ^  »  z1  >  *  •  » 


••  O  « 


iwf  (v  )  -  likvO.  (v  ) 
o  o  —  °  1  0 

5 


0  , 


^104b) 


if  last  equation  is  simpl;  tl  uatioj  of  conti  »ni  ■  1  e  loc 
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re  cor  i  tiering  lo  A.  •  ici  '  ^  <  1  n  o 

in  neglect  the  terns  containing  that  is,  torn  proportional 
to  the  inverse  square  of  t-  v  plasms  frequencies,  (l^f )  then  r-  -duces  to 
the  cimplo  form 
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tho  given  in  (45)  and  (^6)  with  z-  "being  replaced  Tx  , 
that  the  same  relaxation  time  apys-rs  in  tho  ox*  roe  ion  ."or  t'  :  tte- 
nuation  for  both  the  shear  vr  veo  and  th  dilation*?!  •  vo-.  ha  v 

mentioned,  the  expression  (?!)  for  is  just  tl  A  hie’  1 

th-  exprecr  ion  for  tho  sh-  rr  viscosity  of  th  electron  g '  a  i  c  "cul.*  tod 
by  t  inb  rg,  (1950^)  hlso,  our  anal; -sir*  s'  ;r  t-  .  t  th  -  1  ctron- 

phonon  interaction  does  not  load  to  a  bulk  viscosity  of  tha  -loctrcn 
'"’bis  r  ult  is  in  agreement  if'  the  exp^rimen 1  -  vie  an c  a-  of 


gas 


Ilason  and  Bommol,  [4  ]  ,  which  u .  :  the  electro 


bulk  vi  cosi ty,  e  (j) )  „ 


negligible,  if  my, 
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3r  •’T,"*  ■“  p  ATi*JT*;rv  /*  'v't  t*tt\  '  rp  t  »t 
o  J  x  ■  y».  i .»  >  v  i.  X  •  j JL  *  f li 


.  T  r  _  «  \  Y-  >  Mri  tp 

I  -  1  J 


In  section  3.3  end  3.4  v/o  wore  able  to  obtain  ’elutions  by 
taking  only  th<  fir  t  i  in  t1  ■  s  h.  rioal  harmonic  n  ion 

for  f,  ,  (87)*  owevor,  in  V  .  ca  o  of  tl  e  ligher  freon enci  i P  >>  l 

we  rre  no  longer  able  to  terminate  the  'cries  at  a  satisfy  ctory  point, 
'once  n  attempt  a  solution  by  ret'  ini ng  oil  the  terms.  flirthor  in 
all  events,  it  is  necas.  r;  to  obtain  an  expression  only  Cor  J  /e 
or,  equivalently,  Cf  (v0  ) ,  for  both  the  incident  shear  waves  nn<  fco.' 
the  incident  dilations  1  wove  .  (93'  has  been  'written  in  tl  rorm 


(1 


wz: 


..No  rY  > 


r  o 


ik  r. 


p  -l).f  (v  )  +  (p  +  ?)v'L.'7_  fv  ) 
*  p-n  o'  o  pi-  t  o' 


(Pp  -  1) 


O-P  +  3) 


P  =  1,2,3,... 


'  no' 


(llO)  re  pro.  ele  n  infinite  net  of  coupled  linear  algohr  ic  omation. 

as  its  solution  a  set  of  ratios  of  two  infi¬ 
nite  determinants .  It  is  mere  convenient,  however,  to  olve  tl;  fol- 


c  ( v  '  1 

/>'  * 


in  trie  ;  *  s  and 


lov  iny  emotions 

(P  -  1)  ’  p.,^o  )  + 


ia 


P 


P 


w 


i) 


iy,  (p  +  ?)pw(v0) 

Wol 


where 


and 


£E»  (  1  +  ££-£.'  ^  ^>,  *  p  =  1  » 1  ’  * * *  * #  > 


V  7* 


V 


mu 


"v  '  =>  v  b  (v  )  , 

'  n  1  a  Px  o  *  ' 


o  P 


1  h  i*^) 


fill) 
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We  now  une  the  method  of  determinantal  ration  to  write  down 


the  formal  solution  for  this  set  of  linear  equations.  jn  tMf  case 
we  lave  to  deal  wit':  infinite  determinant,  and  so  we  shall  assume  that 
the  ration  converge  to  the  actual  solution  of  +.’  rrcM  Furtl  r, 

we  paint  out  that  n  analytic  solution  over  the  complete  r-  ny*  of  fre¬ 
quencies  has  not  been  obtained  and  .ve  ve  limited  ourselves  to  the 
low  and  high  frequency  regions.  However,  as  di  cussed  Inter,  nur  ri— 
cnl  solution. .  c  n  be  obt.-  ined  fs-r  the  intermediate  frequency  1-  nq  >, 
but  we  -h  ll  not  cornu  rn  ourselv-.  >  with  these  numerical  solutions  in 
this  thesis. 


(ill)  has  the  solution 

P 


r\  / 

•  P  / 


wl  ere 


1 

-lia 

T 

o 

0 


—  >i: 

9 


.o-j 

t3 

0 


0 


7  1 


-iia 

7 


0 


0 


-W  7 

9  J 

1 


0 

0 

c 

-6  is 
11 


D  is  D  with  the  pth  column  repl;  oe  by  roa,/l  f 
f  YZ\  mu 


0 

0 

0 

0 


in  the  first 


rov 


<:  and  aero;  elsewhere  in  the  column.  Tn  this  esse  e  r'-  y  '  ri  I 


J 


( v  ) 


I-  U  A  , 

Q  I  '  o' 

m  v 

o 


where 


// 


-•( 

n  \ 


l  + 

mu 


v-  / 


U  /  1  +  e  <£  '  >  , 

V  mu 


-  / 

I  / 1  JL$V'  . 

mv  V  mu 


o 


1 
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(112)  is  of  exactly  the  same  form  as  (52)  ;nd  thur  th  solu¬ 
tion  can  ho  carried  out  in  the  same  manner.  Thin  re  ;ulto  in 


t  = 


1 . 

im  k  1 


b")  (kV  +  k*  -*V) 

/  V  0  •  • 


and 


I 

c 


IT  mu  ( 1 
o  v 


-  *  ) 


k^b77  +  x 

o 


it 


fence  the  equation  of  motion  of  the  ions  becomes 


8 


/ft 


2)is 


// . 


If  mu  ( l~b  )  (k2-  +  1-  ) 

o  o  7 


(k^b"  +  k  2) 

'  o  > 


(113) 


Again  we  ho ve  precisely  an  analogous  form  for  the  equation:  of  motion 
as  (r>7)  in  the  case  of  an  assumed  relaxation  time. 

// 

To  evaluate  the  attenuation  coefficient  e  crust  evaluate  b 
or  f,/ D,  that  is,  the  ratio  of  two  infinite  determinants,  for  kJ>l<<  1, 
this  ratio  can  be  expressed  in  the  form  of  a  continue*1  fr  ction  which 
can  be  terminated  at  any  point  for  any  desired  approxim  tion.  mVis  is, 
of  course,  equivalent  to  terminating  the  series  expansion,  f  ■  in  section 
3.3,  at  any  desired  term,  and  leads  to  precisely  the  same  results.  W  en 
>  1,  this  method  no  longer  applies*  However,  it  i  found  th*  t  tt 
leading  term  in  b*  can  be  determined  quite  easily.  This  term  is  suffi¬ 
cient  to  give  the  first  order  t  rras  in  the  attenuation  and  T  ■  disper¬ 
sion  relations,  but  does  not  give  the  correction  terms.  7cr  interme¬ 
diate  frequencies,  k f  ~  1,  it  seams  that  one  must  use  r  numericrl  method, 
Because  zp  and  hence  decrease  rapidly  with  increasing  p,  in  pr  ctice. 
for  num  rical  work  v;e  h:  ve  always  only  to  solve  the  r  tio  o  t  o  finite 
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determinants.  For  example,  if  kf  «=  1,  then  YP  =  k  Pt  ^0.01.  ”  nee 

,,  r-  - 

we  expect  that  only  a  10  x  10  matrix  ne  d  to  he  used  fer  this  v"luo  of 
YJ  .  We  shi  11  not  attempt  an;  numeric.  1  solutions  in  t'  i  th<  sis 
hut  shall  asiume  that  the  solution  carries  smoothly  from  t^  o  low  fre¬ 
quency  limit  to  the  high  frequency  limit. 

.  /» 

r"o  evaluate  F,  /V-  end  hence  h  for  kV)>>l,  wc  now  write 


f  =  £3.,/  1  +  c  £  L  , 

Y,z.  V  mu  , 


-7n  *. 


>ia 


5 


0 

0 


0 


iia£ 


~2ia 

5 


0 


0 


v  * 


n 


Dividing  a  ,  out  of  each  rev  we  hove  for  5  ,(v  ) 

r  i  c 
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A  \  V  ^ 

I  ' 


ma  ,/l  +  e  ^ ^ 


kr. 


mu 


0 


0 


0 


0 


m±  0 


-i.  ot3l  ]_ 


-i  * 


jy 


'  3 


•i*a  0 


0 


-1^# 


“i  ^3 


0 


~1«* 


32. 


I_ 


-1  « 


Jy 


where 


f  1 14 


i(<  +!) 


Oi- til  • 
(;i  +  3) 


nd  < 


2i  +  1 


From  the  form  of  the  determinants  v/e  see  that  the  denominator  oen  ho 


expended  in  the  form; 

S  +  ]  Wa/  +  5/i/a 

o  ^ 


•  ••j  +  |  ^2/ /(ai +  4?  ■  /:  3) 


+  •  •  •  j 


4-  •  *  « 


Tie  nurit  rator  can  also  be  expanded  in  t.ie  form. 


1 


<f0+  \  Z&i  +  fa  v  j  +  0  • 0  f 

1  - 


The  ^’g  are  infinite  products  of  prrticvl,  .  <x 1  •  Since 

only  the  leading  torra  in  relatively  simple  to  evaluate  me  U-  rc'1  ;r' 
rather  difficult,  we  shall  restrict  our  calculation  -  to  tv  1  oa<  inf; 
term  alone,  which  allows  us  to  calculate  the  attenuation  coefficient 
but  not  the  correction  terms  in  the  dispersion  rol-tiom  .  Also.  1- 


though  Sc  by  itself  may 


.11  we  need  is  the  exprer,- 


not  remain  finite, 
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aion  for  (v^)  to  remain  finite*  Hence  we  must  evalut  te  2f0  /&  wheri 

L> 

this  term  is  obtained  by  inspection  from  the  determinants  and 


and 


(  (“1  )  (-i^ii)  ... 


~  (-ik3)  (_i  asz)  (**i  5  (_i"^V;  o*. 


:nce 


*0 


5^7)  - 


CO  , 

=  ;Id  f,)Z;} 


30  r  7  .  r 

77"  i*u.-0  *;{}«*,  ti,  -')] 

i  -  I 

t  «.]••<  t'  o  ex’"  ir.  .''hr  1/  c  that  there  is  no 

factor  of  i_  appearing  in  this  term,  "furthermore,  we  see  that  th  .re 
must  be  one  more  oc  appearin.  in  the  denominator  t’>on  in  the  nun  rotor 
as  we  start  from  the  first  row  of  '  for  S  and  the  second  row  o"  -  for 

O  I 

Xp,  To  evaluate  t  is  pr  iuct  31  b '  ct  to 

evaluate  the  product  for  a  finite  number  of  terms  and  f  en  t  ]  the 
limit  as  the  number  of  terms  goes  to  infinity  together,  that  in, 
keeping  the  1)  and  Th  the  same  size,  ^or  r  finite  number  of  terms 

N 


O/N 


JT  \Ul*  ^>0  *(**<■')  Zc 
r/  XFTzT^oJt  *2:Uc  -/)  j 


lAI  1.+.AI 

2N 


N 

.77 

*  -  /  L 


2i  + 


yU  +  3; 


k 


4i  +  1' 


+  1)  (21  +  1 


(4i  -1) 

fc?I  -f) 


(41  + 


1) 


A/  _ 


2N 


IT 

<  -  / 


4i  (i  +  1)  (4i  -1)' 

(21  -1)  (2i  iYJTfttr) 


Uhl 


2 


LUlii) 

2i: 


3  lVff  —  a 

1* a 5 1 71 . . .  T-  -3r(?  -1) 5rr'"  -1 i7( ' 5 ) 
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zn ,  . 

2_Uir  +  l) 

(47  +3)  (27) 


iClL).,  .  .Gi  +  lJ . . 

{(2H-l)l(/4  (2T:+1) 


y  +i) 

r% 

* — s 

%  — , 

_ 1 

(41+3) 

( 2 VJ  )  *■ 

- 

We  now  apply  Stirling's  approximation  in  the  form 


, — —  h? )  -  a/ 

/2W  V  G 


'(’his  gives 


N 


2TI( '.■'••+X)(4r+3) 

-  5fr+x)(-r+i) 

2;;(2?!+i)(4i'+3) 


*/<•  _  -2M 

2  (2<7) ;?  e 

e~2M 


(mi) 


,Aiz 


577Ttt+i  )  (jh+l) 

2  (2r+i)(4K+j) 


N  ->  oo 


2. 


V 


4  (115) 

Therefore  we  have  for  the  leading  term  in  ,T  /o.  from  '  114) 


md  ( 115  )  ^ 

J  , 


=  ’  I  ;  )U  A  f 
mv_ 


rQn  ma ,  (  1  +  e  £  \  1_  & 

mvo  kr(  \  mu  /  a, 


That  is, 
// 


=  K  u  /  1  +  eCZ) 


37T  +  .. 


rii5&) 


mu  y  \  4k 


■>  TT  +  (smaller  term.;); 
4k  ^ 


and  indeed,  h  <<  1. 
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e  have  near  only  to  evaluate  tl  tl  tion  coefficient. 


'■lor  |k*|«|k*|  and  Ik^ « |k’b"|  in  (113) , 


2  k  =  o(  -  (Arm 

X  s 


KV$k' 


■ie  1  -  1 

Vo "  > 


*=  01 


■*  * 

1*1 


S  / 


4k,  f, 

ITT 


4m v 


-  4mv  *>-  Si/) 

Wv}  ' 

This  is  precisely  the  limiting  value  found  in  the  case  in  which  we  as¬ 
sumed  the  existence  of  a  single  relax  tion  time,  (65). 

^or  k^  >  >  |k*b//|  and  k^>>|k^|  ,  the  situation  is  again  the 
same  as  that  found  in  aection  2.3  and  from  (113), 


0/ 


ii 


m 


]  C//^s  ^ 


Re  (1  -  b  ) 


3  M  +  •  •  a 

4k  4 , 


fin) 


Of  particular  interest  is  the  fact  that  only  Tf  now  a  .po  r  in  this 
first  order  expression  for  the  attenuation  coefficient. 


'  .VK 


high  ra  :qu>incy  i  gibs  pl  jil  r  t' 


In  this  case,  we  use  the  same  formalism  pp  developed  in 
section  3«5  and  again  restrict  the  calculations  to  an  expression  for 
the  attenuation  coefficient.  The  same  comments  hold,  hero  as  in  sec¬ 
tion  3*5  with  the  complications  mentioned  in  section  3.4.  We  have 
(l04h)  and  (104a) 


Go(vo)  - 


1  kv^G  f  (vo  ) 

7  t*r 


ana 


(1  +  i*sZ^)C^(vo  )  -  ikr. 


*  0  ,  p  =  0 


P  Vp-/(v0) 


2pirj 


+  (p  +  i)vii „  (v^) 

*7Sv — ; — — s  O  f>+t  v  0/ 

(2p  +  3) 


0(1  +  e£f,p^  ,  p  »  1,2,5, 

imu 


fll8b) 


Again  it  is  more  convenient  to  sot 

A a(v  )  *  v^C-fv.) 

pK  o  '  o  P  o  ’ 

and  solve  tie  infinite  set  of  linear  equations  for  A (v^). 
(I18e)  and  (118b)  result  in 


Ao(vo)  ~  1  ,  (r0) 

3  uT 


0  ,  p  =  0 


(119a) 


and 


-ia, 


p 


4  ( v  ^ 

p  -  <\  o  ' 


+  \  (v  ) 


->  ( p  +  1 )  p  „  ( V  ) 

^WTT)  ^  ' 


ma 


1  +  e£iw  ]£>,  »  p  =  1,2,3, ... . 


k  r/  V  mu 


119b) 


'e  now  have 
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77 


r 


-ia 

0 

0 


c 

3 

1 


-?ia  _ 

3 


0 


0 

-2ia , 

r* 

5 


r\ 

rJ 


0 

0 


-Iis  l 


7 

1 


0 

0 

0 


-41  a  - 

9 


and  D/  is  given  by  replacing  tl  second  column 
cept  for  the  second  row  which  has  the  value 


by 


ma ,  (  1  +  e£jf,\ 
k  Zt  \  mu  J  * 


11  zoro?  e: 


In  E,  ia  *=  kv  and  a  _ 
*  o  —  o  r 


?  ,  as  before.  And 


A 


written 


k  P 

(i  +  i  uj  zyi 


Di  /E. 


To  evaluate  A  f  first  consider  the  ratio  e/t).  •  hicl  can  bo 


2  -  hi'  [  1  +  iLli 


<?zv 


-  / 


ma, 


mu 


la ,  j 

3 


i  j 


v-  ere  '  is  the  minor  of  the  first  element  in  E  and  D.  is  the  minor 

of  the  first  element  in  divided  by  ma,  /  l+e  <f  zt  \ . 

kr(  \  mu  / 


/D.  . 


e  have  now 


to  evaluate 
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t 

T> 

St 

1  'a, 

-i 

0 

0 

0 

«  •  • 

f 

t 

~i  *2t 

1  As 

x'  i 

-i 

0 

0 

•  «  « 

0 

-i 

1 

0 

•  a  • 

•  •  ♦ 

<t  •  • 

•  •  • 

•  •  • 

a  •  9 

1/a , 

-1  ain 

0 

0 

0 

•  99 

0 

l't:x 

-i 

0 

0 

9  9  9 

0 

-1*31 

1/a, 

-i  -*jy 

0 

9  9  9 

•  a  • 

•  •  • 

•  0  • 

•  •  • 

•  «  * 

9  9  9 

where  now 

q/  '  / 

<  ( < 

V/>  - 

i  H*  1 

2i  +  3 

and 

*  i 

?i 

+  1  , 

+  1 

I  he  leading  term  in  a  series  expansion  o 
a  ,  &0  ,  V  where 


will  now  ho 


)  (”i  *t/)  (-1*3*)  (-i ••• 


^  *  (-i*^;  i-i**,; 


and 


once 


O  q  at 

X 


(-i oi23)  (-i«Vy)  -i^yy)  *••• 


[<xz  3sc3i)  ~(*  yy *sv)  7  s*  7^)  *  •  • 


Again  we  sum  for  a  finite  number  of  terms,  using  the  know¬ 
ledge  that  there  are  no  factors  of  i  left  singly  in  this  expression 
and  that  the  product  S  contains  one  moroocthan  does  ^ ,  -inc-  is 
the  term  of  a  determinant  one  degree  larger  than  that  giving  i£.  ”e 
put  this  extra  term  at  the  beginning  of  the  product,  evaluate  Kf^0 


for  a  finite  number  of  term?  and  take  the  limit 
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the  number  of  terms 


goee  to  infinity  together  in  &Q  and  X?  >  as  in  section  3.r). 


o\v 


IT  \ >  .1  v. qza/ 
-/  l<2,  (Z.n)  J 

a/ 

IT 


/V  r/ 


/  -  / 


4i  +  3) 


(4i  -1)  \2i  +  1 


(41!  +  1) 
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(.'»•  +  1)  (4K  +  3) 

y(2K  +:i)ji^ 

5(4!'  +  l) 

- 

r  2/V  i  I 

z 

- 

S5 

A    vJ  v  •  J 

(2H  +  1)  (4N  +  3) 

'  (2!!)J 

3(2::  +  (4:  +  1) 

rx 

77.' (4N  +  J 

_ 

on  applying  .tirling’s  approximation.  In  letting  h  tend  to  infinity, 
wo  have 


'o 


-  77 

z 

Hence , 


■n 


Therefore, 


and 


where 


// 


k  1/  / 1  +  etc, 
ml  mu 


la 
3  ° 


7T 


-  ! 


=  r» 

/ 
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rn 


k  r, 


N  uA  , 

o  ' 


mv 


1  +  e£4/\/^a0  +  77* 

mu  /\  3 


ou  /l  +  efr\  ", 


mu 


o 


7Tk/  /  1  +  2a. 


o 


v 

o 


IT 


+  2kv  \  . 
"IT”  \  i^tr  y 


'120) 


'121> 


'ma) 
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’"e  note  that  (121)  is  in  precisely  the  same  form  as  'llo:  ).  ^he  ana¬ 
lysis  just  follows  that  in  section  2.2.  Neglecting  tl  o  terms  con¬ 
taining  the  plasma  frequency,  the  equation  of  motion  of  the  ion:  for 
k /  >  >  1  is  (analogous  to  that  of  section  2.2) 


3ls 


IT  mu 


/12°) 


?x 


* 


rr 


// 


fsing  (121a)  and  e  uating  imaginary  parts  in  f\  .e  >  ve 


X  =  2k 


in  us  '{e 


V5  k ' 


g 


Therefore , 


//  mv  ^s 
o 


'/5> 


mv 

— ° 

3:/^ 


<X 


'f 


7mt.v  us 
o 

—77- — 


-  / 


7* 


+  mv 


//  mv  ^ 

_ o__ 

^  •? 

OiV  v  7 


•Set?, 


e  note  again  that  (123)  is  just  that,/ (48),  obtained  when 
we  assumed  the  existence  of  a  single  relaxation  time.  Also  the 
term  in  the  velocity  due  to  the  hulk  modulus  of  the  electron  gas  appears 
in  this  calculation  snd 


1  + 


mv 


-° 


\  * 


3’**  J 


4.  :  HI  CU 


rrii  ■  *• 


•  *  «*r  *  *vr  ^  a  j »  rpj  ■  *  ■  ’  ■  j  yy y  y^*** '  f1  r  <  •  f*  *  rm^  wj  ^p*».Tp 


hen  there  ore  impurities  in  a  metal,  the  electron-  are  not 
only  scattered  by  the  thermal  motions  of  the  ions  but  ore  also  scat¬ 
tered  b;  the  impurities.  The  treatment  of  the  collision  integral 
cri  ing  from  the  impurities  is  similar  to  that  for  the  electron-p  onon 
interaction,  (fo) .  Also,  the  rein,  inder  of  the  an  lysis  leading  to  the 
various  expressions  for  the  attenuation  is  suite  analogous  to  th? t  in 
the  previous  sections.  Therefore  vo  shall  not  give  ti  d  ils  of 
this  calculation  and  shell  only  st:  to  the  results.  mhe  only  difference 
is  that  vie  get  the  following  relaxation  times  characteristic  to  e;  ch 


term  in  the  expansion  '07)  , 

,3 


V 


=  (20  + 

27<*. 


■  +  Ip>  f  I  -'M I  [i  -  K  (coon) 

v/!  ere  ft  «  |  K  -  K  ^  (  =  I .  [2{l  -  cor,  A  )J 


'124  ' 

\  1  / 


Air  the  of  cat t 


that  is,  the  angle  between  II  and  K 7 ,  and 


|H&0|*  =  c 

ri(K  -  E  'K(r)d? 

^(20  -i-  1)3 

J 

is  the  absolute  square  of  tie  matrix  element  for  the  transition  of  the 

— *  — * , 

electron  from  the  state  into  the  state  T  0  c  is  th  e  concentration 
of  the  impurity  atoms  and  V(r)  is  the  perturbing  potential  due  to  one 
such  impurity  atom.  In  making  t’  is  calculation  -e  h,  v  assumed  th-  t 
the  electrons  can  be  described  by  plane  wave  func  io  •  it  th.  t  tber ?  is 
only  one  type  of  impurity,  and  that  the  various  impurity  atoms  •  re 
randomly  distributed.  e  remark  that  ox; r  i  n  (124 '  for  p  »  2  ir 
the  same  as  that  obtained  by  ftcinberg,  (l95f3b)  jyj,  in  hi  ca n oula- 
tionr  of  the  viscosity  of  the  electron  gar.. 
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It  is  seen  that  if  |'t(A)|  is  independent  of  the  angle  of 

scattering,  that  is  X,  then  the  various  Z^’s  are  all  equrl  end  just 
*  -2. 

gener<  1,  qti  j  J  depends  on  /l  teen  ^  l ^  m  rro  /^et  some 
idea  of  the  ratios  '//f/  ,  we  have  made  a  calculation  of  |  |  *  by  ta¬ 

king  for  V(r)  a  screened  Coulomb  potential  of  the  type 


V(r) 


~/*T 


~  JiQ 


rl  96) 


where  l/u  in  the  range  of  the  inter-  ction  and  A  in  the  strength  con- 
stant.  :e  find  then, 


m(*) 


)  \  —  C 


1  ''  TT^  i  ^ 

I'  ll  A 


"‘herefore. 


»  1 


TP 


+  i)*  xpr^rp- 

/ 

{l  -  f  cosfl) f  d(cosil)  > 


cos  Try* 


where 


2fo:  1  = 

■a-fi^vK  *  <omvIC 


end 


1  + 

2K' 


(127) 


(1  8) 


1?9) 


'ince^>l,  then  we  see  from  ( 128)  that  Zp  ^  ‘  ,  p  *  1 , 2 , 3 > . . .  > 

V 

for  the  interaction  potential  (126).  The  cose  of  particular  interest 


•  _  i  /  _  t 

IS  Z~  /  4^  o 


'e  may  expect  that  r/ '  will  be  greater  tl  rn  unity  and  will  have 
a  similar  qualitative  behaviour  so  long  as  the  scattering  i  predomi¬ 
nantly  in  the  forward  direction  and  so  long  as  it  is  a  monotoni colly 
decreasing  function  of  fl. 
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pecifieally,  for  p  =  1  ancl  p  =  2,  (l?8)  reduces  t< 


In 


1 


m 


1?0) 


I  •  =  2 


c> 


lnpf  +  1 


^151) 


0  •  Vi  plotted  the  ratio  Z*,Zl  in  figure  4  against  xf,  that  i: 

*  / 

effectively  the  range  of  the  impurity  potv  tial*  "or  ^  ir  large, 


this  is  for  _1  <<1,  »  1.  -hi.-  is  because  tl  e  -c  ‘  '  srin  b  c 

/* 


JL  -  i. 

isotropic  for  this  condition.  On  the  other  hand,  "or  ^  =  0  //  «=  5. 

In  the  case  of  a  polyvalent  impurity  in  a  monovalent  net  1,  the  poten¬ 


tial  is  approximately  of  the  screened  Coulomb  type,  a  ,  ,;hown  b  "ott 

and  Jones,  (l93o)  J^lil  o  hor  example,  the  screenir.  constmt^  is 

-e 

C.55  x  10  cm.  for  such  an  impurity  in  copper.  emerib  ring  that  for 

22. 

copper  the  number  of  electrons  per  unit  volume  i  8.5  x  10  per  c.c. 
then  e  have  =  lo08,  ev  luated  at  the  Fermi  surface.  anc:  ther>  fore 

*0  =  1.9. 

?/ 

e  have  seen  that  for  kl#<<l  the  attenuation  coefficient 
is  proportional  to  Z~2 ,  (10')  and  (108).  If  vie  take  only  th  impuri¬ 
ties  into  account  the  analysis  would  be  preci  .ely  the  rue  nr.  the 
new  expressions  will  be  simply  ( 102 )  and  (108)  .vith  rz  replaced  Iv 
Z'  .  '  Ti on  the  electrons  are  scattered  both  by  impurities  and  by  ther¬ 

mal  phonons  v. e  shall  assume  that  the  effective  relaxation  times  are 
given  by 


P  =:  1  >  ■■  » 3  *  *  • 


frn) 


v 


-p 


r; 
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Actually,  in  our  calculations  (132)  is  not  a.  sepr.ro.te  assumption 
but  follows  from  the  various  assumptions  made  in  sections  3  <r,nd  4. 

For  p  =»  1,  (132)  corresponds  to  Tvlatthiessen* s  rule  for  the  additivity 

reciprocals  of  the 

of  thei relaxation  times  in  electrical  conductivity  due  to  t  impuri¬ 
ties  ••  nd  to  the  thermal  motion  of  the  ions.  "e  nov;  go  on  to  discus' 
the  results  of  this  and  the  preceding  sections  ip  relrtion  to  the  ex¬ 
perimental  data. 
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5o  DISCI!:.;  ION 

To  should  stress  that  our  treatment  of  the  atJenuation  pro¬ 
blem  in  the  previous  sections  is  based  on  two  mein  approximations: 

(1)  The  C'ommerfeld  free  electron  model  of  a  metal  is  applicable  and 

(2)  The  vcrmi  energy  surf.  ee  is  spherical.  In  view  of  t>  c  above  r  .3- 

•U  > 


sumptions  we  cannot  expect  good  quantitative  agreement  with  the  .-hr 
lute  values  of  the  attenuation  coefficient,  particularly  fir.ee  quanti¬ 
tative  measurements  have  be  n  mode  on  lead,  tin,  and  aluminum  only  ,(#J. 
'"ome  measurements  have  also  been  made  on  copper  but  these  } ;  v-  not 
been  corrected  for  the  energy  losses  inherent  in  the  apparatus.  mo 
are  also  limited  in  verifying  all  the  consequences  of  our  work  with 
experimental  data,  for  the  following  reasons:  l)  -’he  calculations 

of  the  relaxation  times  from  the  "Bloch  integral,  at  best,  gives  re¬ 
sults  to  within  an  order  of  magnitude  unity,  as  is  well  known  "rom  si¬ 
milar  calculations  of  the  electrical  conductivity.  (2N  This  difficulty 
could  be  overcome  if  measurements  as  a  function  of  temp  -cture  were 
available  for  both,  the  attenuation  coefficient  and  the  electrical 
conductivity  on  the  same  specimen,  however,  very  few  such  meesur  - 
merit.'  h  ve  been  made.  (3)  Measurements  h ■  ve  not  been  made,  as  *  rule, 
over  a  sufficiently  large  frequency  range  at  a  fixer-  temperature  so  as 
to  cover  both  the  regions  kl/<<  1  and  kl , >  1. 

The  preceding  calculations  indicate  the  following.  rst. 
for  kl  «1  the  relaxation  time  Z  ,  appearing  in  the  expressions  for 
the  attenuation  coefficient,  (102)  and  (100),  in  tie  same  for  both 
incident  dil  tional  and  incident  shear  waves.  ; of erring  to  the  1*  t 
paragraph  of  section  ,  v/e  3ee  that 
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(IV) 


1  1 
Z1 + 


3 


r  1 


r 


where  2T  and  r‘  are  the  relaxation  times  appropriate  for  electrical 
conductivity,  respectively  due  to  tie  electron-phonon  intorcction 
and  to  the  electron-impurity  interaction.  .Moreover,  $  is  a  constant 
of  the  ord< r  of  magnitude  unity*  For  a  typical  impurity  potential 
such  as  the  screened  Coulomb  type  given  in  (126),  6  lies  between  1 
and  3.  '/hen  both  electron-phonon  scattering  and  electron -impurity 
scattering  are  present  the  relaxation  time,  Tc  ,  appropriate  Tor  elec¬ 
trical  conductivity  is  given  by  (132), 


rc 


1  =  1 
r 


4*  /  * 

a 


(X’4) 


r  zr 

Hence  Z^/r  will,  in  general,  lie  between  1  arid  1/  3« 

A  consequence  of  the  fact  -t hat  Tu  is  the  same  for  both  the 
dilational  and  the  shear  waves  is  that  the  ratio  ,  (102\  and  , 
(108),  is  given  by 


'135) 


This  ratio  seems  to  be  verified  by  experiment,  [4 J  ,  for  single  lead 

in  the  norm  a  i  state 

crystals^  t  2  K  and  for  *  2(  .5MC/sec. ,  oi  «  O.rv  r  6;m. 

r 

V  =  2.39-10  om 'sec.  For  at  =  10.lI.IC/sec. ,  *  *  0.070  nepers /cm.  and  m 

^  /  3  1  3  $ 

3  2 

1.27xlO’rcm*/sec.  fence  =  0.67  and  Jj/VA  /^-A  -  0.73.  Measure 

4  Vv5  /  \*rt) 

ments  of  the  attenuation  of  both  dilational  and  shear  waves  on  other 
substances  are  not  yet  available. 

Oecondly,  a  method  of  comparing  our  theoretical  results  on 
the  relaxation  times  with  experiment  would  be  as  folio  s.  First,  we 
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fit  the  theory  to  the  experimental  data  at  the  low  frequency  limit  by 
multiplying  our  expression  (102)  with  an  arbitrary  cor  t  nt.  b.  pi- 


-  / 


versus  k  jP 


;uro  5  shows  the  plot  of  j\*.,  {J) /*rj  lira  ^  [kr)  'ur^ 

*•*  c<r->oo  ^ 

for  measurements  on  very  pure  tin  crystals,  ^4*].  In  Pine  ro's 
c;  leu  lot  ion,  vvhere  T  is  identified  us  2^  ,  the-  eight in  '  ctor  is 

about  1  ",  resulting  in  the  curve  indicated.  In  our  c  ue  • -n  w- 

*• 

•  r  «  Z^/3,  the  weighting  factor  is  about  3/2  nd  he  c 
high  frequency  limiting  v.  lue  will  be  just  three  times  t ■  t  n?  pip- 
pard .  Z~c  ,  ’  ,  i  termined  by  electrical  conductivity 

measurements.  Cur  limiting  value  of  the  dil  tionel  attenuation  *t 
bight  frequencies,  is  much  closer  to  the  experimental  value  tv  ■  n  r'ip- 
pnrd’a.  If  we  take  the  impurities  into  account,  t}  on  it  i  cl <:  r 
from  (133)  and  134'  that  lies  between  Z“c  nd  £~c  • 

improve  the  agreement  with  experiment. 


A  further  consequence  of  (133)  and  (134)  is  that  the  r;  tio 
Ta  t'Tc  is  temperature  dependent.  fence  the  weigl  tin  factor,  b, 
will  depend  upon  the  tamper  ture  and  so  also  the  limitin  value  •  t 
high  frequencies  of  the  r;  tio  plotted  in  figure  5»  On  t’  e  b  sis  of 
our  work  thi  limiting  value  will  decrease  ;itl  decreasing  t- m  .na¬ 
ture . 

"e  note  that  shear  waves  behave  somewhat  differently  than 
do  the  dilations  1  waves.  If  we  fit  the  attenuation  for  kfj  <<  1,  as 
for  the  case  of  dilational  waves,  the  v<*.  lue  for  the  attenuation  t 
^  1  will  again  be  about  2*5  time,  the  value  obt  ined  b  sott.inv 

M. 

her euroments  by  llson,  (l?39)  1.14]  ,  give  vr lues  of  thi  'lighting 
factor,  b,  of  3  for  tin  md  4*5  for  aluminum. 
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Za  *  lc  ,  throughout.  ?or  k^>>l  and  » 1 ,  however,  the  attenu¬ 


ation  approaches  the  limiting  vnlue 


m  which  will  he  about  throe 


times  that  expected  for  TA  *  T£  .  It  ■  to  tt.'  in 

the  condition  <^^>>1  requires  much  hi  ;her  frequencies  tl  h  v  y  t 

b  on  obt-  iried  experimentally « 

throughout,  we  have  con  idered  proces  es  hi ch  con'  ribute 

both  to  the  attenuation  and  to  the  electric.'.  1  rc  i,  torico.  Th  r-r  i  •, 

however,  one  process,  namely,  the  electron-electron  colli  ion-  •■•■hie1 

males  no  contribution  to  the  electrical  resistance  inoe  the  total 

electron  momentum  is  conserved  in  these  collisions .  the  other 

hand,  these  collisions  will  contribute  to  attenuation  ju  t  ■  the 

atomic  or  molecular  collisions  contribute  to  viscosity,  jn  -■  oners  1 

in  good  metals,  because  of  the  strong  screening,  the  electron- electron 

inter.' ction  will  not  be  of  much  importance.  If  they  are  import- nt, 

as  may  be  the  case  below  1°K,  then  v/e  will  have 

1=3  s  1  fUG) 

£ a  £"/  Z~,  t-ee 

er<  is  tl  •  rel  x  tion  time  ]  propr  t  1  ctron-electron 

-Z 

collisions  and  is  expected  to  be  proportional  to  for  rmi- 
Dirac  gas.  1/k  will  then  have  the  form 


I  “ 

-x 


•»  '/n  Z. 


n 


Myy ) 


whereas,  the  electrical  resistivity,  ,  would  be 

f  »  A’  i  M7.A) 

Hence,  by  plotting  l/«i  end  />  f  r  the  same  pecimen  oi  •  r  y  obtain 


inform  tion  cone  rnin 


90 


Although,  as  already  mentioned,  the  electron-electron 
collisions  are  not  likely  to  be  of  import.- nee  in  metals,  t '  ey  may  be 
important  in  some  semi -metals  arid  semi-conductors,  ^he  att- rust ion 
measur  ements  may  th-en  be  used  to  obtain  inform*  tion  core  min  the 
relaxation  times  for  the  electron-electron  colli  ion;  , 


Value  of  kft  . 

(1)  Measured  value  for  very  pure  tin  at  60  megacycles. 

(2)  Pippard‘3  theoretical  results. 


(3)  Our  theoretical  results. 


Figure 


-  / 

)/<~r  lim  fcXy,  (ter)  /^]  1 


I  • 


Plot  cf 


versus  kJ 


PART  II 


ULTRASONIC  ATTENUATION  IN 
POLYCRYSTALS 
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6.  INTRODUCTION. 

As  we  have  mentioned  in  the  general  introduction.  Mason  and 
McSkimin,  (1947)  [l5_|  »  (1948)  [  16  j  ,  have  shov/n  that  one  of  the  me¬ 
chanisms  through  which  attenuation  arises  in  polycrystals  is  the  scat¬ 
tering  of  sound  waves  by  the  grains.  The  scattering  arises  from  the 
fact  that  the  effective  elastic  constants  relevant  for  the  propagation 
of  these  waves  differ  from  grain  to  grain  because  of  the  different 
orientations  of  the  crystallographic  axes  of  each  grain  with  respect 
to  a  coordinate  system  fixed  in  space.  In  the  frequency  range  for 
which  the  wavelength  is  large  compared  to  the  grain  size  Mason  and 
McSkimin  find  that  the  amplitude  attenuation  coefficient,  *< .,  ,  for  both 
dilational  and  shear  waves  in  polycrystalline  aluminum  and  magnesium, 
may  be  written  as 

-  V  +  V*  (!) 

The  first  term  is  attributed  to  an  elastic  hysteresis,  for  which  the 
precise  mechanism  is  not  yet  known,  and  the  second  term  to  scattering. 

In  this  part  of  the  thesis  we  shall  deal  with  a  calculation  of  I)  « 

When  a  plane  dilational  wave  is  incident  on  a  region  of  mat¬ 
ter  whose  elastic  constants  or  density  differ  from  those  of  the  sur¬ 
rounding  bulk  matter,  both  dilational  and  shear  secondary  waves  are, 
in  general,  produced.  If  the  sheer  modulus  is  zero,  as  for  a  fluid, 
only  the  former  are  produced.  Mason  and  McSkimin  realized  that  in  a 
solid  both  types  of  scattered  waves  are  produced.  However  in  analy¬ 
sing  their  experimental  results,  they  used  a  formula  calculated  by 
Rayleigh  which  refers  to  the  scattering  of  sound  waves  in  a  fluid  medium. 
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The  value  of  obtained  is  several  times  less  than  the  observed  value 

✓ 

for  polycrystalline  aluminum. 

In  the  case  of  incident  shear  waves,  Mason  and  T'c^kimin  use 
a  formula  which  is  written  down  in  analogy  to  Rayleigh’s  formula  for 
the  scattering  of  light  due  to  inhomogeneities  in  the  refractive  index 
of  the  propagating  medium.  The  variations  in  the  refractive  index  of 
the  medium  arise  out  of  density  variations.  t.Vq  note  that  the  scat¬ 
tering  of  light  is  due  solely  to  the  density  variations,  as  evidenced 

by  the  angular  distribution  of  the  scattered  intensity.  On  the  other 

sound 

hand,  in  solids,  thej^scattering  is  primarily  due  to  the  variations  in 
the  elastic  constants.  However,  by  identifying  ,  the  density,  in 
Reyleigh's  formula  with  yu  ,  the  shear  elastic  constant,  they  obtain  a 
value  for  which  is  about  the  same  as  the  experimental  value  for 
aluminum. 

Since  the  formulae  used  by  Mason  and  McSkimin  are  such  as 

to  have  incident  plane  dilational  waves  scattered  only  into  dilaticnal 

waves  and  incident  shear  waves  scattered  only  into  shear  waves,  Rhatia, 

(1959)  [17] ,  has  improved  on  their  calculation.  His  work  is  bared  on 

the  following  assumptions:  (a)  The  grain  size  is  small  compared  to 

the  wavelength  of  the  incident  wave  and  that  all  the  grains  are  nearly 

effective 

the  seme  size,  (b)  The  changes  in  the^elastic  constants  and  in  the 
density  of  the  medium  from  gr-.in  to  grain  are  small,  so  that  they  may 
be  treated  as  a  perturbation  and  multiple  scattering  can  be  neglected 
and  (c)  The  grains  are  randomly  oriented,  so  that  the  bulk  material 
may  be  treated  as  elastically  isotropic  and  characterized  by  Lame’s 
two  elastic  constants,^  and  yt^  .  A  further  simplification  was  made, 
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this  being  that  the  effective  elastic  properties  of  each  grain  are  des¬ 
cribed  by  just  two  Lame’s  elastic  constants,  A  and^*,  where  \  -A0«  K 
and  /<.  -/*0<</u£  His  solution  consists  now  of  solving  for  the  scattered 

energy  from  a  single  grain,  summing  over  the  scattered  energy  from  all 

* 

the  grains  in  an  element  of  volume,  and  hence  obtaining  a  value  for  the 
attenuation  coefficient.  The  main  result  of  his  work  is  that,  for  both 
incident  dilational  and  shear  waves,  more  than  90y  of  the  scattered 
energy  is  oarried  away  in  the  secondary  shear  waves.  The  attenuation 
coefficients  obtained  now  are  about  twice  too  large  for  the  incident 
dilational  waves  and  only  slightly  larger  for  the  incident  shear  waves 
than  the  experimental  values.  The  theoretical  ratio  of  the  two  atte- 

reasotiubiy 

nuation  coefficients  is  about  six,  agreeing^well  with  the  experimen '.al 
value  of  15.5  for  aluminum,  but  not  with  that  of  unity  for  magnesium, 

[15}. 

Actually,  an  elastically  tonisotropic  grain  has  several  inde¬ 


pendent  elastic  constants,  their  number  depending  upon  the  crystalline 
symmetry  of  the  grain.  As  a  consequence,  there  ere,  in  general,  twenty- 
one,  (not  necessarily  all  independent),  non-vanishing  elsstic  constants, 
C!  .,  referred  to  a  set  of  orthogonal  axes,  the  laboratory  coordinate 
system,  not  coincident  with  the  crystal  axes  of  the  grain.  mhis  part 
of  the  thesis  takes  the  existence  of  these  elastic  constants  into 
account  for  the  attenuation  problem.  Our  assumptions  are  just  those 
of  Bh&tia’s  with  (c)  modified  so  that  all  the  elastic  constants  of  the 
individual  grains  are  retained.  The  solution  follows  the  same  pattern 


as  his 


wre  also  assume  that  the  clastic  anisotropy  is  small. 

In  section  7  we  formulate  the  general  equations  of  motion 
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for  an  isotropic  bulk  material  containing  a  single  grain  and  solve 
for  the  scattering  amplitude.  These  solutions  are  valid  for  polycry¬ 
stals  containing  grains  of  any  crystal  symmetries.  In  section  8  we 
connect  the  expressions  for  the  energy  scattered  by  a  single  grain  to 
the  attenuation  coefficient  of  the  sound  wave.  The  resulting 
expressions  for  attenuation  involves  certain  averages  and  mean  devia¬ 
tion  averages.  These  averages  are  computed  in  Appendix  I  for  crystals 
of  orthorhombic  symmetry.  Crystals  of  higher  symmetries  are  now 
simply  special  cases.  Finally,  in  section  9  we  discuss  our  results 

in  relation  to  previous  theoretical  work  and  the  experimental  observa- 

* 

t ions. 


*The  results  contained  in  this  part  of  the  thesis  have  appeared  under 
the  title  Scattering  of  High  Freauency  Sound  Waves  in  Polycrystalline 
Materials”  in  Jour.  Acous.  Soc.  Am.  ^1,  1140  (1959). 


7.  FORMULATION  AND  SOLUTION  FOR  THE  SCATTERED  LAVE  FROM  A  SINGLE 
INHOMOG EN  El  TV . 

The  equations  of  notion  for  a  solid,  in  the  absence  of  body 


forces  are: 


a* 

St 

+ 

dx 

£U  + 
ay 

ilf. 

/-A 

•  It* 

sr 

1 16  * 
ax 

*JLz  + 
dy 

an  ^ 

^  d  w 
'  at* 

ss 

U.S  + 

y  t  ^  4. 

d  X 

dy 

where  yo 

is 

the  density,  u, 

v,  and  w  are  the  Cartesian  components  of 

the  displacement,  and  the  T^’s  are  the  components  of  the  stress  tensor. 
For  small  displacements  it  is  a  consequence*  of  Hooke’s  Lav;  that  tho 
stresses  are  proportional  to  the  strains  and  hence  v;e  have 


where  we  use  the  convention  that  repeated  indices  are  summed.  The  C. . 

ere  the  elastic  constants,  and  the  S.'s  are  the  stresses  defined  bv 

J 


S 


dX 


a  v 

a  y 


d  w 
d  z 


V 


d  w 

ay 


Bv 

dz 


S  - 
j 


~  ^ w 
dx 


3u 
a  z 


s 


c>  V 

a  x 


+  d  u 

ay 


(4) 


It  is  a  consequence  of  the  fact  that  the  internal  energy  is  a  perfect 


differential  that  C.  .  «*  C.. 

id  ji 


and  hence  there  are  at  mo^t  21  indepen¬ 


dent  elastic  constants ,  decreasing  in  number  (to  two  for  an  isotropic 
crystal),  as  the  crystal  symmetry  incretses. 
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7.1  TEE  EQUATION"  OF  MOTI  N  FOB  THE  BULK  MATERIAL. 


In  general,  one  denotes  the  elastic  constants  by  C.  when 
the  axes  of  a  fixed  coordinate  system  coincide  with  the  crystallogra¬ 
phic  axes  of  the  single  crystal.  Usually,  this  situation  is  not  ac¬ 
complished  and  the  single  crystal  axes  are  oriented  with  respect  to 
the  fixed  axes.  In  this  event,  the  equations  of  motion  in  the  fixed 

coordinate  system  .are  given  by  (2)  with  the  C . . ’ s  replace  by  C! ..  the 

^  J  J 

transformed  elastic  constants0  The  C . . ’ s  depend  onlv  on  the  C  ’  s 

*  ij 

and  on  the  angles  between  the  fixed  axes  and  the  crystallographic  axes. 
These  angles  may  be  described  in  terms  of  the  Euler  angles  for  the  ro¬ 
tation  of  the  crystal  axes  from  coincidence  with  the  fixed  axes  to  their 
actual  position. 

The  bulk  material  is  composed  of  small  grains  arbitrarily 
(or  randomly)  oriented.  Since  we  are  considering  frequencies  such 
that  the  wavelengths  of  the  ultrasonic  waves  are  large  compared  to  the 
dimensions  of  the  individual  grains,  we  assume  that  the  zero  order 
elastic  constants  for  the  bulk  material  are  the  average  values  of  the 

t 

C  ’  s  over  all  random  orientations.  This  results  in  the  bulk  material 
J 

being  described  by 


Co  =  C',\)  “  (Cn)  ”  (p's?) 


1(  C„  +  C*2  +  C  ) 

3 


-  2P 

15 


a  -  <c;y>  "  <CA>  -  <ci«> 

=  l(  C  wy  +  C  -  -  +  C.J  +  lr_  )  ^5) 

_ 3 

^  The  transformed  elastic  constants,  the  Euler  angles,  and  methods  of 
averaging  are  explained  fully  in  Appendix  I  and  the  results  only  shall 
be  stated  here  and  in  the  following  sections. 
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-  1(C  ,  4  C  +  C  )  +  IP  , 

5  "  15 

and  ^\Cij/“  0  *  for  1  =  1»2>3  811(1  J  =  4,5,6> 


where  P 

=  +  +  °7s')  ~  (  + 

C/j’  + 

C«) 

~  2(  C*,  + 

tv  ♦ 

That  is,  we  consider  the  host  material  as  an  isotropic  medium  with 

the  elastic  constants  given  by  (5)  and  the  equations  of  motion  given 

by  (2),  'Ye  denote  the  solutions  to  these  equations  of  motion  by  uQ, 

v  ,  and  w  . 
o'  o 

We  shall  assume  that  u,v,  w  and  u, ,  vq,  are  proportional 
to  a  simple  harmonic  function  of  time,  exp(iu/t),  where  «/-  =  2J7f  is 
the  angular  frequency  of  the  incident  wave  and  u,  v,  and  w  are  the 
solutions  to  the  generalized  equations  of  motion.  The  time  dependence 
will  be  suppressed  except  where  noted. 

'Ye  may  now  define  the  propagation  constants,  (=  2  7T/A), 
k  and  h  for  the  dilational  and  shear  waves,  respectively,  as 

k  as  =a  US'  ,  h  =*  ^  ^  j  [  6) 

C  V4  s<0  V% 

where  V,  and  V,.  are  assumed  to  be  the  zero  order  velocities  of  the 
dilational  and  shear  waves  in  the  bulk  material. 


•*The  values  V,  =*  and  V  =  f/<o  are  expected  to  be  in  agreement 

with  observed  values  to  within  only  about  10%  because  of  the  extreme 
simplications  made  here. 


7o2.  THE  EQUATIONS  0?  MOTION  PON  THE  INHOMOG HN12T T Y . 

’7e  consider  now  a  single  small  grain  arbitrarily  oriented  in 
this  bulk  material  to  be  impinged  upon  by  an  incident  wave.  Tn  the  inho- 
mogeneity  the  equations  of  motion  are  given  by  (4)  where  we  must,  of 
course,  use  the  transformed  elastic  constants.  These  we  write  in  the 
form 


■  cij  +  scij> 


(7a) 


where  the  C^’s  are  the  average  values  given  by  (5)  and  where  we  shall 
assume  that  the  S C.  ,’s  are  small  quantities,  that  is,  they  may  be  trea- 

**“  J 

ted  as  a  perturbation.  ’7e  also  assume  that  the  solutions  of  the  equa¬ 
tions  of  motion  can  be  written 


u 


u  +  &u,  v 
o 


V  4-  Sv.  W  »  w  +  &W 

o  *  o 


'7b) 


Hence  we  may  also  write 


j  =  h  +  SV 


(70) 


restituting  (5).  (7a,U  (7b),  and  (7c)  into  (2),  retaining 
only  linear  terms  in  all  the  small  quantities,  because  of  our  assump¬ 
tion  of  small  anisotropy,  that  is,  in  all  the  s,  and  making  use  of 
the  fact  that  \ic  ,  vc  ,  and  wc  are  solutions  to  the  average  equations 
of  motion,  we  obtain 


/? 


or  &u  + 


a  \  c  s  s 
-  0  1 

+  A/s s,  +  gs  ) 

ax 

J 

— 

e)Z  i  J 

"  "  [fy>  + 

2>y  t 


/«£  '6 


(8a) 
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4  J_j>„SSsl  +  _i  j  CoSS^  +  ^CfS,  +  SSj)] 


dx 

^y 

a  z  l 

«  -  8/>  cs v  +  H 

y  °  y  - 

9 

(8t>) 

A,"*6*  +  -2.USsj- 

dx 

!  +  J.  +  -± 

ay  L  -1  hz 

c  ss, 

0  J 

+  A0(SS,  + 

S  SJ)  |  *  -  S/O  us1  wq 

+  H  / 

TV 

4J  . J 

(8c) 

where 


H 


i  (SC,  ,s°  )  +  _a  (Sc6j3°  )  +  JL  (Sc,,s®  )  , 


9  x 


ir  j 


dy 


az 


5J  j 


T? 


—  ^C6jSj  ^  +  ~  -is-i  )  +  jL  )  i 


a  x 


ay 


2  j 


3z 


y-j  j 


(8d) 


H 


_2  (sc  s®  )  +  j_  (Sc  ,s°  )  +  ±  (sc  s®  )  , 

ax  ^  ay  J  J  dz  J 


summed  over  j .  We  note  that  the  left  hand  side  contains  only  the  per¬ 
turbed  portion  of  the  wave  and  the  average  values  of  the  elastic  con¬ 
stants  and  that  the  right  hand  side  is  determined  completely  by  the 
initial  wave  and  the  inhomogeneity* 

Pi "ferentiating  each  of  (3)  by  x,  y,  and  z,  respectively, 
and  adding  we  have,  after  making  use  of  (4)  and  (6),  and  setting 


s  m  ui  v  •+  wk  ^ 


and 


H 


Hi  f  Kj  +  H  k, 

irJL  7.— ’ 


•  &  a)  +  C  y^(V»S's)  =  -  [^i\7*(^8o)  +  \/-TT  J 

or  \/\\/  •  &  s )  +  k  (  v*5«f)  =  -  1^ \/« \$fks2  +  n  j 

C0  L  ° 

-  -  V'f/ c0  > 


(9) 


10? 


where  V  is  the  usual  Laplace  operator  and  9  is  the  usual  divergence. 

Differentiating  (8a)  with  respect  to  y  and  (8b)  with  respect 
to  x  and  subtracting  and  continuing  this  process  we  obtain  the  three 
equations $ 


\7  (curl&s  )  +  h^  (curias  )  =  -  _1_  curl  j  8/?  ur2z  + 

L  X  ° 

=  -  curl  * 


do) 


The  solution  to  (9)  now  represents  the  dilational  component 
of  the  scattered  wave  which  is  produced  when  the  incident  -wave  passes 
through  the  inhoroogeneity ,  and  the  solution  to  (10)  represents  the  shear 
component  to  the  scattered,  -wave.  Here,  we  shall  consider  two  sets 
of  initial  conditions, the  first  is  a  plane  dilational  wave  propagated 
in  the  negative  x  direction.  That  is, 


s  =  u  i  =  Aexp  (  ikx  )  i  .  (11a) 

o  or“  — 

The  second  is  the  case  in  which  we  have  a  plane  shear  wave  propagating 
in  the  negative  x  direction  polarized  in  the  y  direction.  That  is. 


c 

o 


•lib) 


*  =  Aexp  (  ihx  )  j.  -• 


7.3  SOLUTIONS  FOR  TITO  SCATTERED  WAVE. 

The  particular  solution  to  (9)  which  represents  a  suherical 
outgoing  wave  is  given  by 


=  +  _ 1_  f  e  '  Vc.yf(xo9yo,zo)&ro  , 


4^C  J  r 
o 


in  which 

z 


r"  =  (  X  -  Xo  +  (  y  -  yo  }i  +  (  z  -  Zo  ^ 


JL  i 

dx 

o 


3  j  +  <3k 

dy  dz 

^  0  o 


(15b) 


\ 

d  z:  =  dx  dy  dz  , 
o  o  *  o  o  ’ 


fl3c) 


end  the  integration  is  over  the  whole  volume  of  the  solid.  In  this  cc.se 
the  integration  goes  only  over  the  volume  of  the  inhomogeneity .  n*  nee 
the  S’ s  and  their  derivatives  are  zero  everywhere  except  in  the  inhomo- 
geneity  and  on  its  surface. 

Integrating  (12)  once  by  parts  and  applying  the  above  boundary 
conditions,  that  the  S  ’ s  and  their  derivatives  are  zero  outside  the  in- 
homogeneity,  and  further  using  the  fact  that 


JL 

3x 

o 


t 


we  obtain  for  (12), 


L 


(14a) 


*  The  integration  of  (12)  follows  that  given  by  Rayleigh,  Theory  of  Sound, 
(  Dover  Publications ,  New  York,  1945  )»  Vol.  2,  p.  149. 
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In  exactly  the  same  way,  the  solution  to  (l6)  is  given  by 


curlSs  =  curl 


L  ^ycJ  r 


e~‘  r  V"(xo*yo«zo)dr» 


Cl4b) 


In  (14a)  and  (14b)  we  have  integrals  of  the  type 


J  »  e  (x  ,y  ,z  )  Sz>  (x  ,y  ,z  )dr 

l  —  o  o’  o  o'  /  v  o’  c*  o' 


o 


'15a) 


and  JT  a 

2 


''  e~4/<rr  3  (5c.  .  (x  ,y  ,z  )s°  (x  ,y  ,z  )\dc 

-  \  1J  ^  o^o'  oy  J  v  o’“o’  oyj  o 


(15b) 


to  evaluate.  Integrating  once  by  parts  ana  applying  the  boundary 
conditions  above,  then  it  reduces  to  an  integral  of  the  type  J  .  That 
is,  it  becomes 


A  "  C  A  /  —  J  V  O  /  \ 

d  e  §C..(x  y  z  )S.(x  ,y  ,z  )dr  . 

- —  /  —  ijv  o^o*  o'  y  o’^o’  o'  o 

ax,j  r  °  ° 


(15c) 


Remembering  that  the  volume  of  integration  is  small  compared 

f*'  ^  . 

to  a  wavelength  of  the  ultrasonio  wave,  then  9  x  S,<x  ,y  z 

~  y  o  ©  © 

will  be  a  slowly  varying  function  over  the  region  of  integration  and 
hence  we  may  replace  it  by  the  value  at  the  center  of  the  inhomogeneity. 
In  particular,  if  we  take  the  origin  of  the  coordinate  system  at  the 
center  of  the  inhoraogeneity ,  then  we  may  set 


x 

o 


z  =  0  * 

o 


and  hence 


J 

t 


»  £  4 ><>  uTl  s^(0,0,0)  />dro 
r  ^ 

*  £  ‘  ^  ^^^(OjOjO)  /U?  T  * 


(17a) 
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where  now  r  is  given  by  (lja)  and  (l6) ,  T  is  the  volume  of  the  inhomo¬ 
geneity  end  £yo  is  the  average  value  of  over  the  inhomogeneitv.  In 
solids,  it  turns  out  that  y?  can  be  considered  to  be  constant  over  the 
volume  of  the  inhomogeneity  and  hence  iy>  =  Cy?  if  we  neglect  boundary 
effects,  which,  in  fact,  will  be  negligible a 
In  the  same  manner, 


d  x 


Sj (0,0,0)  /  ^  i - (Xq *0 o » 


dx 


zlkr  G°(0,0,0)T£C 


ij 


(17b) 


where  C  is  the  average  of  <$C  .  over  the  volume  of  the  inhomogeneity . 
i  j  i  j 

Our  model  gives  SC  to  he  constant  over  the  volume  of  the  inhomoge- 

ij 

neity  except  at  its  surface.  For  all  practical  purposes,  these  surface 
effects  are  negligible  in  AC.  and  we  may  set  AC  =  6C  . .  '’’his  is 
done  throughout  the  rest  of  our  calculations.  Further,  we  note  that 
7,  given  by  (l3a)  and  (16),  the  outward  vector  from  the  cent  r  of 
the  inhomogeneity  along  which  the  scattered  wave  propagates. 

Using  (4),  (8d),  (9),  (14a).  (14b),  (I7f),  and  (17b)  to¬ 
gether  with  ecch  of  (lla)  and  (lib)  we  have  the  solutions  for  the  in- 
cident  dilational  waves  and  the  incident  sheer  waves j 
(i)  The  incident  plane  dilational  wave  prepresented  by  (lla); 


V-&S  -  V9  J£L  I  ^A^Ri  +  ik^/ 


4  77Co 


(13a) 


and 


curl  S  s  =  curl  aT  j  w'2  A^>  qi  +  ik  '/¥t 

AW 


(13b) 
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where 


-  (  ac„  *2  +  4c/6  11  +  AC/caK)  i  +  (  /SC  3_1 


H 

it 

4  lO 

c'Ar 

» 

— * 

wnd  F 

1 

AC,,  11 

( 

d  X 

A 

C,  dR 

'  ay 

+  AC/y. 

Q.  = 


-  i'hr 

®  * 


dy 


15 


dz)  ( 


'6 


dx 


%S1  *  A0,* 21)1  +  (  CiO  21  +  AC  an  +  ac  ai)k 

9y  az)  (  dx  dy  adz) 


lSc) 


(I8d) 


end  ^  is  simply  H/  with  R  replaced  by  Q. 

(ii)  The  incident  plane  shear  wave  represented  by  (lib); 


\/'S  s 

=  (7.  AT  +  ihiq  ]  , 

4  "  c0 

(19a) 

and  curlSs 

=  curl  AT  i  o*r2H/o  Q  j  +  ihA(  j  > 

4  /7>0 

(191)) 

*  —  y 

where  the  quantities  are  all  as  above  except  that  H  and  are  H, 

and  IV  with  AC  .  replaced  by  AO... 

'  'J 

At  distances  far  from  the  inhomogeneity  we  can  neglect  A/r 
compared  to  unity  ana  the  operations 


_a  (A‘kr\ 

dx  (r  ) 


-<  Ar 

-  il<:Le 


r 


etc . 


7e  define  the 


direction  co  ines  cf  the  vector 


as 


,  M  -  £  ,  N  « 

r  r 

That  is  L,  M,  and  IT  are  just  the  direction  cosines  of  the  outward  going 
vector  from  the  inhomogeneity . 

The  solutions  (18a)  and  (l8b)  can  now  be  rewritten  as 


V‘6 


8  -  !'  +  k^/°" 

4/7’ C  " 

0 
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(21a) 


and  curl  S  s  =  curl  ATQ  i 

4.^  >- 


+  khn: 


(21b) 


where 


11°  =  (  Lac,,  +  Mac,s  +  nac/y  )i  +  (  Lac ,  + 

HOC,*  +  NAC/y.  )i  +  (  Lac,f  +  Mac^  +  VM,3  )k  ,  (21e) 

I 

Similarly,  the  equations  (19a)  and  19b)  can  now  be  rewritten 


as 


and 


v-$- 


4  7TG  L  ' 
o 


“SI  . 


(22a) 


curl S s 


curl  ATQ 

4  ^ 


+  h*  ^  j  , 


(22b) 


"jb  “3  -*  -* 

where  IH  and  Ha  are  related  in  the  same  way  as  IT  is  to  TH  . 

to  16 

Finally,  to  the  first  order,  given  by  (20),  the  solutions  of 
Sq  at  distances  far  from  the  inhomogeneity  oan  be  put  in  the  form 


Ss  -  S,  +  .  (23) 

The  solutiorefor  the  incident  dilational  wave  ere  obtained  from  '21a) 
and  21b)  resulting  in 


and 

=  \TR  j  (ur*  ISfi  i  + 

4TTC  L  " 

0 

,  2rO)  — 
k  II,  )  *  n 

4_ 

n  j, 

-j 

s' 

J 

=  -  _ATJ  j  {u-3hp i 

+  khH,  )  x 

n  j  X  n  j 

(24a) 


'24b) 


where  n  is  the  outward  going  unit  vector, 


n  »  Li  +  M£  +  ?Tc  „ 


(25) 


The  solution*  for  the  incident  shear  wave  are  obtained  from  (22a)  and 
(22b)  resulting  in 


4ffc 


ATE  "[  +  WiH^) 


n 


AmQ  j  {“S~Z n  j  X  n  - 

4^ 


(26a) 


In  the  above  expressions  we  see  that  represents  a  spherical  outgoing 
dilational  wave  propagating  in  the  direction  of  the  radius  vector  at 
a  velocity  and  with  displacement  along  the  radius  vector.  Similar¬ 
ly,  S,  represents  a  spherical  outgoing  shear  wave,  propagating  with  velo¬ 
city  Vs  along  the  radius  vector.  Its  displacement  is  perpendicular  to  the 
radius  vector  and  its  polarisation  is  a  function  of  the  angle  of  scat¬ 


tering. 


8. 


TKE  ATT  CITATION  COEFFICIENT. 
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TTIK 


SC  ATT  BRING 


CROSS-SECTION  FOR  A  SINGLE  IKTIOFOO  LIT  ST  TT  . 


We  define  the  scattering  cr oss^aection,  ,  as  the  energy 

scattered  per  unit  time  into  a  unit  solid  angle,  centred  r.t  the  scat¬ 
tering  center,  divided  by  the  incident  energy  passing  unit  aroa  perpen¬ 
dicular  to  the  direction  of  propagation  of  the  sound  wave  per  unit  time. 
First,  the  total  energy,  due  to  the  incident  wave ,  (lla),  in  an  element 


of  volume  dV  is 


dE  =  IdV  =  1  LvA  |  />c  d?  ;> 

r> 

c 


where  I  is  the  intensity  of  the  sound  wave,  that  is, 


(27) 


X 


2  2 

2 


/° ? 


( 28) 


If  we  have  a  volume  Y  with  unit  area  perpendicular  to  the  direction 
of  propagation  of  the  sound  wave  and  with  length  Y  ,  then 


A 


=  1 


ft 


UT 


l\2  y 

-*  v 


(29) 


represents  the  incident  energy  passing  unit  area  per  unit  time. 

Similarly,  the  energy  passing  unit  area,  perpendicular  to 
the  direction  of  propagation  of  the  sound  wave,  per  unit  time  due  to 
the  incident  sound  wave,  (lib),  is 


w  J  =  1  ■  (30) 

2 

Also,  the  energy  passing  an  element  of  area,  dA.  on  the  sur¬ 
face  of  a  sphere  with  the  inhomogeneity  at  its  oenter  due  to  the  spheri¬ 
cal  outgoing  scattered  dil&tional  and  shear  waves  is, 
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no 


5  w 


>  V  +  %VS  )  dA  , 


2 


(31) 


where  now  and  pj  are  the  maximum  values  of  S^,  and  Ss  ,  respectively, 


2 

and  dA  =  r  sineded^. 


Fence,  by  definition, 


^-(e(jzS)dA  =  _1_  (  kV  .+  S2V  )dA  , 

.  A*VS 


(32) 


and 


cr(e9/)dn 


_i_  (  +  Sh’s  )«»< 

A%  *  ' 


(33) 


where  da  =  sinsded^  and  we  use  (24a)  end  (24b)  for  and  $T.  in  ^3?), 
and  (2 6a)  and  (26b)  for  S*  and  S5  in  (33)  • 

The  total  cross-section,  <7~  ,  is  now  obtained  by  integrating 
o~(e over  all  angles  of  scattering,  This  integration  is  quite  straight¬ 
forward  and  expressions  (32)  and  (33)  need  only  be  written  out  fully. 

The  result  for  (52)  is 


7 


T 

477-0, 


1  rT  1  +  1)2‘['|  (AC  ,  f  +  (AC^f  + 


3 


^  ^ 

(AC  3)  )  +  £(  («  ,)  +  (AC  r )  +  (AC„)  )  +  2 (  (00  )X 

)  3(  ) 


(AC^)  +  (AC„  )  (&Pl3)  +  (AC/<2)  (AC,3) 


3( 


+  V,  JZ 

Tt  W 


£Jty>)‘ 


_2k  h 
15 


2,-i  K  (AC„  )*  +  (AC  f  +  (A  Z  3f  )  +  3(  (AC  f  +  (AC  ) 


/-T 


+ 


)  ( 

(ac„  )  (ac/a)  +  (ac7/)  (ac^)  +  (ac^)  (ac^)  > 


7  Z 


/  / 


Rewriting,  this  expression  c;n  be  put  in  the  form 


_ri-“ 

4^/ 


(1(4.)  +  jfB)  +  Ma.  (  10V)  +2 

(3  V*  )  }/(3  7 


OI ^ 


k  h  A 


5^ 


Ill 


4y?T*f* 

r 


1(  1  +  ?/V/\3  )  (Q*f 


L3( 


>  /S 


I,  (  E,  +  2/1- 


5Co2( 


r 


4<Ay  £o 

Y*V  ‘ 

/  •<? 


1  +  1  /V. 

2\w 


+  2 

5c2 


j/ 


A 

■A.  , 


+  1/Ys\ 

?(  Y  j 

VV  / 


^  *  ) 
'  \ 


(34) 


In  the  Sr.. me  way,  we  obtain 


<7~ 

5 


„,2  ,T 

Y/ 


4/IC«  Y 

£>  O 

3  2.  *  i 

i  m  a  /»  ^ 


5  5 


+  T 


2  r  ? 


4>/ . 


_2  ^  (zy>)  ■+•  2h  A 


i£H"  £o  f  Sfi?  (  1  +  1/Vi\3  )  (M  +  _2  (  A  .  +  1/ViY  3  )  i 

v  v  j- 3C/ (  2ly)  ;  I*  sc/  (  skj  % 


(35) 


In  (54)  rnd  (55), 


I(  1  (-  i2)  +  (“ci j)  '  *  ')  (aCi6)  +(acij) 

K*)')  -I(  (dCif)  (ac12)  +  (AC±/ )  (ACjj)  +  (acu)  (ac^)  ) 
J  2( 


and 


(36) 


Bi  -  (  (*c1()  +  (ACla)  +  (ac.j)  )  +  i(  Uci6)  +  (ACi5.y 

(  ;  5 


+  (ACivJ 


+  |[  kCi/)  ^°i2)  +  kCi/>  k°iJ>  +  fe'lj)  J, 


(37) 


We  note  that  (34),  (35),  (36),  and  (37)  apply  to  all  crystal 


symmetries.  However,  we  have  limited  the  following  calculations  to 
crystals  of  orthorhombic  symmetry. 


8.2  THE  ATTEITIJATION  COEFFICIENT  FOR  A  POLY  CRYSTAL. 

If  our  only  scattering  mechanism  is  that  of  varying  orientr- 
tion  of  the  grains,  then  we  may  calculate  a  value  for  the  attenuation 
coefficient  from  the  preceding  model.  To  do  this  we  need  tc  make  the 
following  assumptions:  (l)  All  the  inhpmogeneities  are  very  nearly 
the  same  size,  T.  (2)  Each  grain  scatters  independently  of  the  others 
(assumption  of  random  orientations  of  the  grains),  and  (3)  The  varia¬ 
tions  in  the  density  from  grain  to  grain  are  uncorrelated.  The  corres¬ 
ponding  assumption  about  the  elastic  constants  is  contained  in  (2). 
These  assumptions  enable  us  to  sum  over  the  scattered  energies  from 
the  individual  grains. 

Remembering  the  definition  of  the  attenuation  coefficient, 
it  is  the  reciprocal  of  the  distance  in  which  the  intensity  of  the 
sound  wave  drops  to  l/e  of  .its  value.  Hence  it  is  defined  by 

o<  dx  =  -  dl  ,  (38) 

I 

or 

I  =  IQexp(-o<x)  ’ 

where  I  is  the  intensity  at  the  reference  plane  x  =  0,  I  is  the  inten¬ 
sity  at  x,  and  dl  is  the  change  in  the  intensity  in  the  distance  dx. 
From  (29)  we  see  that  dl/l  is  equal  to  the  energy  scattered  out  of  the 
volume  Sax  per  unit  time  divided  by  the  energy  passing  the  area  S  per 
unit  time,  where  3  is  an  area  taken  perpendicular  to  the  direction  of 
propagation  of  the  sound  wave.  Hence, 

___  <  /«  -  Sdx  <C7~  > 

all  grains'-  0  3  T  S 

in  Sdx 


dl 

I 


) 
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(39) 


-  -  l<cr->dx 

T 

where  T  is  the  volume  of  a  single  grain,  and  (cr)  is  the  average  value 
of  cr  over  all  grain  orientations. 

Comparing  (30 )  and  (39)  then 
oC  =>  Co  3  , 


Trom  (34)  sxid  (35)  vre  have 


and 


(41) 


The  averaging  of  and  Th  over  all  grain  orientations  is 

carried  out  in  Appendix  T  find  the  results  are  stated  below. 


8P"  +  _8_  (  4^  +  2b  +  3c  +  d  )  } 
W)  135 


=  JLZ2  +  JL  (  12a  +  b  +  4c  -  2d  )  , 
150  90 


(42) 


where 


114 


a  ss 

+ 

c,-; 

+ 

)*  - 

3  (C^Cjy  + 

<v°« 

+  ) 

b  S3 

(°//  + 

c« 

+  cJi 

Y  - 

3(C/(C^  + 

c/,c^ 

+  c,, c^) 

c  * 

2  + 

c,? 

+  ^23 

Y  - 

3(C/.2  C/j’  + 

+  c,3c.»> 

d  S3 

c,,  (c 

lZ  + 

2CW 

)  +  C^(C/4 

+ 

-  2CV?) 

+  +  C,.5  "  2cW  -  (43) 

Expressions  for  C^,^,  \&l)  >  and  appropriate  for 

crystals  of  higher  symmetry  than  orthorhombic  may  now  be  obtained 
as  special  cases  of  (42)  and  (43  )•  For  example,  we  have  to  put: 

(i)  For  isotropic  crystals; 


I 

// 

“  C  S3 

*=  C  =  ^  +  2yC<  ) 

** 

■  °JV  ” 

C« 

=*  * 

*  °/3  * 

C<? 

% 

II 

(44a) 

(ii)  For  cubic  crystals: 


c//  - 

"zz 

C?3  ? 

c#f ■ 

<Vr  ■ 

c/2  - 

c/j  “ 

^2  3  # 

(44b) 

For 

hexagonal  crystals; 

C  = 
// 

c«  * 

C/J  “  ”  Csr  > 

= 

K  c„ 
2 

“  C,i  )  >  C3J  *  C/Z  - 

(44c) 

^  ,  (40),  and  }  (4l),  represent  the  intensity  attenuation 

0  O 
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coefficients.  Since  the  intensity  is  proportional  to  the  square  of 
the  amplitude,  the  amplitude  attenuation  coefficient  is  just  one-half 
that  for  the  intensity.  Hence 


x,0 

It 

*1% 

(nepers/cm. /cycle  )  , 

(45a) 

d  -  = 

(nepers/cm./cycle^)  . 

(45b) 

9.  DISCUSSION 


The  terms  in  the  attenuation  coefficients,  (40)  and  (41), 

arising  from  the  variations  in  the  density  are  exactly  those  obtained 

by  Bh&tia,  (1959)  i IT  j  •  He  has  discussed  these  terras  fully.  Doth 

* 

the  variation  in  density  from  grain  to  grain  and  the  variation  in 
density  in  an  individual  grain  due  to  thermal  agitation  are  negligible 
in  solids. 

Neglecting  these  density  terms,  then 


(46a) 


(46b) 


The  first  term  in  the  square  brackets  arises  from  the  scattered  shear 
wave  and  the  second  term  arises  from  the  scattered  dilational  wave. 
Since  \  ~  l/2  for  most  solids  and  (A ^  ,  generally the 

second  term  is  usually  only  a  few  percent  of  the  first.  This  accounts 
for  the  discrepency  in  Mason  and  McSkimin1 8 ,  [l5j  ,  [l6j  ,  theoretical 
value  for  the  dilational  waves  as  they  take  only  the  second  term  into 
account.  (Of  course,  they  realized  the  limitations  of  their  calcula¬ 
tions.  ) 

In  crystals  of  cubic  symmetry, 

<"A/>  =  <b6>  =  ^  “  22^*  O3,  /  “ 

6  25.  *  50  75 

where  p  =  C  -  C/Z  -  2Cyr.  Hence  , 

We  note  that  our  D's  correspond  to  the  B's  appearing  in  [16  !  and 

[17]  o 
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(47) 


Ills)3  V  =  6,  f ?  £  *  '  1  +  g/Vs^  ]  , 

4 \yJ  '  12nw  L  3\T/  j  - 


Since  Vs /v,  ~  l/2,  then 


T z s  6  .  ^48) 

(48)  is  just  the  same  ratio  as  that  obtained  by  Rhatia,  [  17  I  .  how¬ 
ever,  we  see  now  that  (47)  is  valid  only  for  crystals  of  cubic  symme-  1 
try.  for  crystals  of  other  symmetry  D  / depends  on  the  values 
of  the  elastic  constants  themselves  through  the  ^  At*^  *  s  and  the  s, *  s. 

The  agreement  between  the  theoretical  and  the  experimental  values  of 
D25/D  ^  ,  evaluated  in  Table  2,  is  reasonable  for  aluminum  but  is  rather 
poor  for  the  magnesium  sample. 

’"e  mention  that  the  measurements  have  been  made  on  two  alumi¬ 
num  samples  having  different  average  grain  diameters.  The  dependence 
of  the  attenuation  coefficient  on  grain  size  is  in  accordance  with 
(46)  to  within  five  percent.  The  dependence  of  attenuation  on  grain 
size  is  the  same  in  all  three  theories;  Mason  and  KcSkimin’s,  i  16  1  , 
Bhatia’s,  i  17  ,  and  ours.  '.Ye  note  that  and  I>2i  depend  quite  cri¬ 

tically  on  the  average  grain  diameter,  as  the  third  power. 

Finally,  Table  2  contains  the  quantitative  values  for  / T 

and  Vis/?  for  the  aluminum  and  the  magnesium  samples.  As  we  have  al¬ 
ready  mentioned,  Mason  and  McSkimin  obtain  values  much  too  small  for 
I  ^  and  close  to  the  experimental  value  for  .  Bhatia,  on  the  other 
hand,  now  obtains  values  in  quite  good  agreement  with  the  experimental 
values  for  both  and  D^s  .  Our  values  are  several  times  too  large 

for  both  D  and  ,  thus  making  the  agreement  poorer  with  experiment 


o 
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In  the  case  of  magnesium,  the  discrepency  between  experiment  and  the 

x 

various  theories  is  much  larger.  It  appears  to  us  to  be  somewhat  pre¬ 
mature  to  comment  upon  these  discrepancies  in  view  of  the  following. 
’Irst,  there  are  very  few  experimental  data  available  -  only  one  set  of 
data  for  magnesium.  Secondly,  as  we  have  seen,  the  attenuation  depends 
sensitively  on  the  grain  size  which  is  not  always  easy  to  measure  accu¬ 
rately.  There  are  also  the  various  simplifying  assumptions  that  we 
have  made  in  our  calculations ,  in  particular,  the  random  orientations 
of  the  grains  and  small  elastic  anisotropy,  which  may  not  be  complete¬ 
ly  applicable  to  the  experimental  data  available.  Renee,  more  experi¬ 
mental  measurements  would  be  desirable. 


**  Mason  and  Mcokimin,  (1947)  Ll5 J  >  quote  a  theoretical  value  for  D2S 
for  magnesium  which  is  in  very  good  agreement  with  their  experimental 
value.  However,  they  obtained  this  agreement  by  using  a  vslue  for 


<AC^)  ~  4  x  10'^,  where,  in  fact,  this  ratio  is  about  4  x  10 


At, 
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Mason  & 
McSkimin 

Bhatia 

Ours 

Observed 

Aluminum 

.  'L*f 

D^/T  x  10 

0.096 

1.29 

4.51 

0.61 

D„/T  X  10** 

13*0 

10.6 

38.2 

8*2 

DiS/D^ 

135 

8*3 

8J; 

13.5 

Magnesium 

/ 

D^/T  x  10 

0.17 

1.06 

9*0 

•  •  • 

DJS/T  x  102f 

20.2 

12.5 

42.8 

0.23 

D2JS  /*W 

119 

11.8 

U.7 

1 

TABLE  II 

Values  of  D^/T,  D^/T,  and  for 

aluminum  and  magnesium* 


APPENDIX  I.  METHOD  OP  AVERAGING  OYER  RANDOM  GRAIN  ORIENT  ATT  OT7S 


The  elastic  constants  transform,  under  coordinate  trrnr- 
formations,  as  the  components  of  a  fourth  rank  tensor  and  for  the 
purposes  of  transformations  it  is  convenient  to  relabel  them. 
Plence,  in  this  section  we  write; 
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II 
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» 
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C33/J 
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Z13,3 
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C*<  = 
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'23 >2  * 

cfjr  ■ 
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~/3i3 
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'/3/Z 

> 

C£6  “ 

C/2/Z  * 

Rirther,  we 

i  note  that 

V 

'  C  .  . 
Ji 

and 

that 

c 

ijkl 

=  C.  ...  a 
ijlk 

C  . . .  -  =  c 
Jikl 

jilk  ' 

A.  2 

The  transformed  elastic  constants  are  given  by 

» 

C ....  *  a,  a .  a,  a.  C  )  A.  3 

ijkl  im  jn  ko  Ip  mnop 

where  again  the  repeated  indices  mean  summing  over  that  index  and  where 

the  a,  ’  s  define  the  rotation  of  the  crystal  from  an  orientation  in 

im 

which  its  crystallographic  axes  coincide  with  the  laboratory  coordinate 
system  to  its  actual  orientation  in  the  solid.  This  rotation  may  be 
defined  in  terms  of  Euler  angles  ^  */3i  ^rid  and  the  rotation  is  car¬ 
ried  out  in  the  three  stages:  (l)  rotate  about  the  initial  z-axis  by 
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the  angle  c*  in  a.  clockwise  direction,  (2)  rotate  abort  the  new  x-axia 
by  the  single^?  in  the  clockwise  direction  end,  (5)  rotat'  abort  the  new 
z-axis  by  the  angle  #  .  These  rotations  are  illustrated  in  figure  6, 
This  transformation  is  given  by; 


n 

=  A 

/  X  \ 

S3 

/a'' 

a  /z 

A!3  \ 

\ 

l 

x  \ 

y’ 

y 

4/ 

a  22- 

a.?J  1 

! 

1 

y 

\  z'  I 

\  f 

u.  j 

\aj, 

Q.JZ 

aS3  J 

\zl 

A.  4» 


with 

A«=  /cos«:  cos  X  -  sin  ex 


;CO?«CUSO  —  3in<*  COS^/f 

oos^xsinb  +  sino(cos^?  cosjf 
sin  c*  sin^ 


sin#'  —  sin<*  cos  S'  -  coso<  cos^  siny  sin^?  sinS^ 
-  sine*  sinks'  +  cos<xcos^?  cor.#'  -  sinxycosJT 
cos  si iyj 


cos^ 

A.  4b 


where  0  <  ^  2  7f  ,  0  <  /3  ^TT  ^  and  0  ^  ^  2  ^  .  A,  4c 

To  average  the  quantities  in  A, 3  over  all  angles  of  rotation 
we  need  tc  v/rite 


/  o  \ 
x  ijkl/ 


1 

QJT3 


J 


r  hjusin/*^fio<dar  ’ 


A. 5 


where  the  integrations  go  over  the  whole  range  of  the  angles  involved 
and  where  the  element  sin^?  d^7  d ot  is  the  solid  angle  produced  by  the 
first  two  rotations  and  the  is  the  element  of  angle  over  the  se¬ 
cond  rotation. 

In  the  case  of  the  average  over  the  product  of  four  elements 
from  the  rotation  matrix,  A,  we  have  the  following  cases  with  thoir 


values ; 
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(a)  All  four  a’s  are  the  same; 


/a*  \ 
\  ira  7 


=  1/  5  - 


A.  6a 


(b)  Two  a's  ere  the  same  and  two  are  the  same,  excluding  (a).  There 
are  nov/  two  possibilities; 

(i)  Both  from  the  same  rov;  or  column; 


<x 


*  „  z  \ 


/ 


7  °  ‘  a  .z  ^ 


im  in  '  N  im  jm 

(ii)  Both  from  different  rows  and  columns; 


=  1/15  * 


A.  6bi 


<aimajn  > 


=  2/l5  »  A.6bii 

(c)  Three  of  the  a's  are  the  same  and  one  is  different.  In  all  such 


cases  all  the  averages  are  zero. 


=  0  <■  A .  6c 

(d)  Three  different  a’s  with  one  souared.  All  averages  are  zero. 


y  2  x  /  2 

<a.  a.  a.  >=  a . 

x  1m  in  io7  im  jm 


in^lan  ^  ~  ^  aima  jn^o  7 


=  0  *  A , 6d 

(e)  All  four  a’s  are  different  then  we  have  only  one  case  of  a  non¬ 
zero  average  and  that  is  when  the  a’c  in  the  product  form  a  minor  of 
some  element  in  the  original  transformation  matrix.  A,  and  hers 


(a.  a.  a.  a.  y 
x  in  im  jn  jm  7 


- 1/30 .. 


A.  6e 


To  obtain  the  desired  average  values  it  is  necessary  to  ex¬ 
pand  A. 5  into  its  full  form  and  substitute  the  above  averages.  When 


this  is  done,  we  have 
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C 
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<C/,„  '  =  <°2  ZtX> 
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15 


A.  7a 


■where 


p 

“  (C///,  ♦  0 2222  +  ?J3333  ) 

^Z//ZZ  + 

+  22  33 ' 
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+  2323 ' 
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'  //  ZZ  33  '  K  /Z 
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-  2(G^  +  <lr.r  +  )  ;  A. 7b 
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l(C/^ 
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+  C /-?/./  +  ^  +  -P 

15 
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l(c  + 

T  ** 

v +  c«  >  +  #  , 

15 

A.  7c 

\  - 

<c/,« : 

>  =  ^(-//33/  ~  V  ^ 2233  's 

S3 

/  /  ZZ 

3 

+  G  //y  j  +  2  33  )  +  JLP 

15 

8 

Kc/z  + 

+  )  +  — Lp  - 

A.7d 

5  15 


These  relations  are  valid  for  crystals  with  21  elastic  constants  and 
can  he  reduced  to  apply  to  crystals  of  higher  symmetry  by  use  of  equa¬ 
tions  (44a)  to  (44c)  in  section  8*2.  '-Ye  note  that  Cq  *■  Xc>+  2 m0 
as  is  the  case  for  isotropic  crystals  and  that  only  nine  of  the  origi¬ 
nal  21  elastic  constants  appear  at  this  points 


These  expressions  reduce 
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to  tho3e  calculated  by  Mason  and  McSkimin  for  crystals  of  cubic  sym¬ 
metry.  All  other  averages  turn  out  to  be  zero. 

To  evaluate  the  average  values  appropriate  for  section  8.2 
we  need  to  evaluate  expressions  of  the  type  <^AC^)  (AC^  )/•  In  our 

model  AC.  «=  6  C .  and  we  have  set  ✓ 

lm  im 

SC.  =  c!  -  c° 

ira  im  im 

from  (9) .  Hence y 


I  f 


_  ~  _T  ’  o  _o  © 

=  C.  C.  -  C.  C.  -  C.  C.  +  C.  C®  , 

im  in  im  in  im  in  im  in  9 


and  therefore  } 


<(ACim)  (ACin)> 


x  »  »  \ 

<c.  c.  > 

N  im  m  x 


p<3  pG 

im  m 


A. 8 


»  i 

Thus  we  are  left  to  evaluate  the  averages  of  expressions  like  C, 

im  in 

or  in  tensor  notation, 


»  t 

C4  m  -.C  ,  ,  =  a.  a.  a,  a  a,  .a  a,,C 

ijkl  abed  im  jn  ko  Ip  ae  bf  eg  dh  r 


j  _ 


A.  9 


where  all  the  repeated  indices  are  summed  over.  Thus  we  have  to  cal¬ 
culate  the  average  values  of  all  possible  products  of  eight  a’s.  All 
the  possible  values  are  tabulated.  In  the  following  repeated  indices 
are  not  summed  and  i  /  j  /  k  /  1  and  m  /  n  ^  o  f  p. 

(a)  8  alike; 


(  a  6  \ 

Vaim  ' 


(b)  7  flike  and  1  different; 


A.  10a 


All  zero 


A.  10b 
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(c) 

6  alike  and  2  alike; 

<  a  a  y  =  ( a .  a  .  ; 

^  im  in  7  x  im  7 

ss 

1/63  , 

*a.  aO  y 

N  im  jn7 

ss 

4/63  - 

A.  10c 

(d) 

6  alike  and  2  different; 

All  zero  , 

A.IOd 

(e) 

5  alike  and  3  alike; 

All  zero  , 

A.IOe 

(f) 

5  alike,  2  alike,  and  1  different; 

All  zero  - 

A.IOf 

(s) 

5  alike  and  3  different; 

All  zero,  except 

<  a.  a  a.  a.  \ 

\  in  in  jm  jn  / 

C5 

-  1.  1  . 

2  ^3 

A.10g 

(h) 

4  alike  and  4  alike; 

\a.  a,  >  =»  \a.  a.  ) 

ira  in  x  im  ora  y 

3 

1 

105  ' 

/  +  \ 

\  a,  a  .  > 

imjn  x 

S 

2 

35 

A.IOh 

(i) 

4  alike,  2  alike  and  2  alike; 

y  V  ^  ^  \  /  *V  .2  J?  \ 

^  “im^ jm^km  /  =  ^aiir:aineip  ^ 

ae 

1 

315  ' 

^aimajmajn  > 

85 

, 

315 

y  #  i  Z  \ 

C  a ,  a  .  a .  / 

x  im  jm  jn 

S3 

— ^ j 

315 

S3 

315  ' 
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/  2.  2.  V'  v 

\Aimajm®lai  S 

,  ¥■  Z  2  \ 

\a ima jn^kp  2 

0)  4  alike,  3  alike,  and  1  different; 

All  zero  . 

(k)  4  alike,  2  alike  and  2  different; 
All  zero  . 

(l)  4  alike  and  4  differnnt; 

All  are  zero  except , 


(ni)  3  alike,  3  alike,  and  2  alike; 

All  zero  * 

(n*  3  alike,  3  alike,  and  2  different; 

,33  v 

Ca.  a.  a.  a.  ) 
x  im  jrn  in  jn  y 

<  a .  a.  a .  a  .  > 

x  im  jm  in  jn  x 

and  all  others  are  zero* 

(o)  3  alike,  2  alike,  and  3  different; 

,3  2  \ 

(a,  a,  a.  a,  a,  > 
x  im  in  jm  jn  ioi  x 

,  z  ?  \ 

\  &,  a ,  a  .  a ,  a  .  ; 

Kp  im  jm  in  jn  x 

/  J  2  v 

(a.  a.  a.  a.  a.  ) 

^  im  in  ip  jm  jp  ^ 


6 

H5 


y 


i L. 

315 


_  1 


2  315 


1.  JA 

2  315 


2  105 

I*  _JL  , 

4  105 


-  I*  _I_ 
4  7 


>15 


-  1.  13 

4  315 


-1.1 

I  315 


A .  lOi 


A.IOj 


A.  10k 


A. 101 


A.  10m 


A.IOn 
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^  aimaina jnakmakn  ^ 


-  1.  1 


4  315 


(p)  3  alike,  2  alike,  2  alike,  a  d  1  different; 

All  zero  * 

(p)  3  alike  and  5  different? 

^aimaina jma j  paknakp  ^ 


b  _L, 

4  315 


Any  one  is  arbitrarily  cubed  and  all  the  others  are  zero, 
(r)  2  alike,  2  alike,  2  alike,  and  2  alike; 


<a?  a?  a  ?  a  .z  \ 
x  lm  jn  jp  / 


<! 


2  2  Z  2  \ 

a,  a.  a.  a.  > 

ira  jm  in  jn  x 


af  a,4  \ 

\  im  in^ai  icp  / 

at1  af  cl-2  \ 
\  ira  in  jp  j£p  / 

/  2  Z  ?  0  \ 

A  a.  a.  ar  a,*  f 
v  im  in  jn  ko  x 


315 

"  I*  -L.  , 

2  315 

“  b  JL  , 

2  315 


315 

b  JL  , 

2  315 


(s)  2  alike,  2  alike,  2  alike,  and  2  different; 
All  zero  , 

(t)  2  alike,  2  alike,  and  4  different; 


c\ aimaina jroa jnakpa jp  / 

/  1  Z  \ 

\a.  a.  a.  a.  &.  a,  > 
•  im  in  ip  jn  jp  icp  x 


All  the  others  are  zero. 


(u)  2  alike  and  6  different 

^-aimajmakmahjiakpainajp  / 


\ 


-  __1_ 

315  ' 

-  1.  1 

4  315 


315 


A.IOo 


A.IOp 


A.  lOq 


A.  lOr 


A.  10s 


A.  lOt 


A.  lOu 


129 


(v)  8  different  3 

All  zero  . 


A.  lOv 


To  evaluate  the  average  of  A, 9  over  all  crientations  one  must 
write  the  expression  out  fully  for  each  case  of  interest  and  use  equa¬ 
tions  A. 10a  to  A.IOv  to  write  down  the  results.  For  our  purposes  we 
need  the  following  -  which  are  written  cut  for  crystals  of  orthorhom- 
hie  symmetry  and  where  we  have  set 


T, 


„  1  nz 

znn  +  c 


1221  +  ZJJ3  3  5 


T 


n,  2  n  CL 

Z2323  +  Zl3i3 


+  Z,2iZ  ’ 


+  Z1/ZZ  +  Zt'/33  +  CZ233  ’ 


'////  ^2212 


+  C////  Z3333  +  C 2222  Z3  3  33  3 


V  =  C/ / / /  °i/22  +  C///>  Z  1/33  +  Z2222  Z//2Z 


+  Z2Z22  Z2233  +  Z3  333Z//3J  +  3333  °2Z3S  9 


T 


I/ll  "22  33  +  Z222  2  Z I  /  33  +  C 33  33  Z  H  21  9 


T 


'/  /  /  /  Z 23 23  +  Z 2  2  2  2  L/3  13  +  Z3333  C'/2/Z  * 


TB  =  °////  Z!3t3  +  G////  Z/2/2  +  Z2222  C/z/ 2 


+  c. 


2222  ^2323  +  3333^/3/3  +  ^33  33^2323 


+  C, 


T?  ~  Z//22  Z//33  +  '  H  22  2233  +  °// 33  °22 33  * 


T 


IO 


Z//Z2  Z 2323  +  C  /  /  ^  -2  Z/3  /3  +  Zf  /  33  Z2323 


+  Z/!33  Z/2/Z  +  ZZl33  ]/2/2  +  Z2l33Z/3l3  9 


T 


11 


+  +  Zz233  Z2323  * 


Z  1/22  Z/2/Z  T  ^  1/33  ^/3/3 


™/2  ~~  '  23  23Z 13  /  3  +  Z  2323  ^/2/2  +  ~!?)3  y,J2/Z  'j> 


A. 11 


<»«..> 


iiii 


1  \  35T  +  43T.  *  12T .  +  6T  +  20T  +  4T 

JjJ  L  '  z  J  **  ^  t 


8T7  +  40To  +  8T-  +  16t  -r  49T . .  + 

/  b  7  /  O  / 1 


32TJ 


<ciii  > 


3L-  f  3  Tf  +  24T2  +  39Tj  +  3T^  +  12T^_  +  6t 


315 


+  CT7  -  12?^  4-  30Trt  -  12T_  -  36T 


T  ,<rn 


<5 


/o 


'// 


or 


/2 


/c’  o’  \ 
\iijj  iikk  / 


1 

315 


T.  4-  OT.  4-  15Ta  4-  3T^+  10m  +  14T  . 

'  -c  J  7-  b  £ 


-  8T7  -  4Tg  4-  52T9  -28T_  -  I2T/y  4-  l6T,„]> 


// 


/-Z  J 


1 

315 


3T,  -  ftp,  -  3?^  -  OT,  ♦  3Tff 

2  2 


-  6t7  4-  6Tp  +  ot^  4-  CfY/0  -  6t;/  4-  nr 


iii jCi jkk  ^ 


315 


8  ■  ^  t  v/-  /<?  '-'-L//~r  /z  \} 

I  [  T  -  6T  4-  2T  -  IT  4-  4T  -  7T 
R  l  '  7  6  7 


<cuu  > 


+  4T7  -  Tg  -  2Tf-  T/a  +  OT,,  +  4T/4 


-i-  |)t,  +  12T2  +  3T,  -  3?y  -  1f  -  2T, 
315 


4T7  -  2Tg  -T,  -  2T„  +  12T„  -  4T„  j 


>  ■ 


1 

315 


T7  4-  12T^  4-  12T^  4-  OT^  4-  T^_  -  4T^ 


+  4T7  -  4T^  -  HTf  4-  8T/<?  -  12T/7  -  8T/X 


> 


1 

315  L 


3T  4-  27T  +  6T  +  3T  -  6t  4-  OT 

’  s/  V  J* 


4-  oT>  4-  6t„  4-  OT.  -  6T  -  6t  4-  !.6T 
?  8  9  /o  "  n 


<C  > 


315 


T  4-  9T,  4-  3T.  4-  OT  -  2T  4-  2T 

9  z  J  *r  o 


-  2T„  -  T„  -  2T  .  4-  2T  +  CT  ,  -  %\  „  |  ;> 


7  8 


/o 


// 


/  z 
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A.  12a 


A.  12b 


A.  12c 


A.  12d 


A.  12g 


A.12f 


A.  12* 


A.  12h 


A.  12i 


m 


i 

315 


(~  5T.  -  24T,  +  12?  4-  6t  4  26?  4-  lOT^, 

l  1  *  3  S  6 


+  2T^  -  2Tg  +  20?  y  +  16T/o  +  12T„  -  16t,2  ,  A.12J 

— 

Applying  A. 8  to  (36)  and  (37)  ©f  section  8.1  then 

<Ai>-l[{0+<.0+  <C^>j-  {(C“,)2  +  (c»/  +  (C»3fj 


+[fci*3+  <'Cir3  +  +  (ctr^  * 

-i[{<</4>+  <C1,0+<CUC^  “  Kcu<,cl<tcU 

and  is  Siven  by  the  3arae  factors  but  with  different  numerical 

coefficients . 

Using  (A*l),  (a.7*)  to  (A.7d),  and  (A. 12a)  to  (A.12j)  then 

f.O  +  O +  O]  -  [<ox +  o2  +(£>* 


1575  L 


128?  4-  312T?  +  193T  -  94T.  +  24T^  -  116t, 

•  J  r  O 


-  72T7  -  32Tg  -  164T7  +  72T/o  -  8T//  -  176T,  1  3 


[ <C > +  <0 +  <OJ  -  [(ftf  +  O2  +(02J 


1 

515 


11?  4-  36?  4-  18?  -  6T  -  T  -  8T 

/  Z  J  *T~  5  Q 


4T7  -  8Tg  -  13?  ^  +  4T/<3  +  12T  /  -  16t/7 


[<p,’,0  +<c)  o  +  <c;2c;_,> 


Q°  £°  +  p©  pi®  ,  p©  p®  ' 

//  /2  //  /3  /Z  /J 


1575  L 


6?y  -  156T  -  29?  -  23T  +  58T  -  82T 


+  36T7  +  16T^  4-  12T?  -  56T/o+  4TX/  +  88T^  ;  , 


s 


132 


Henoe  3 


(a  )  =  _1_  r  952T  +  33*T  +  l65f.T,  -  532t  -  i6*t 

1575  L  3  *  3  J  3  ^  3  s 


and 


<B/, 


// 


-  511T^  -  56t„  -  56t^  -  1252t  +  56t  +  56t  -  168t 

j  6  7  Q  y  9  / o  n  iz 


205^,  +  448T  +  299§T?  -  149^,,  +  56t,_ 
3  '  Z  3  3  3  9  S 


,  A. 13* 


1575 


-  224Tg  -  74fr7  -  74|^0 


Similarly, 


242fr9  +  Ufto  +  74fr„  -  224t/z  , 


]- 

A.  13b 


<\C7 6  +  +  “  [^76^  +  +  ^*5)  J 


1 

315  L 


11T  +  36T  +  18T  -  6?  -  T  -  8T  -  4T  -  8T„ 

l  Z  3  V-  3T  6  7  8 


-  13^  +  4T/<?  +  12Ty/  -  16T  ^ 


[<<">  +  <0  +  <c*>]  - 


«$,>"  +  (c;6^  ♦  (°;£) 


O  \2-  v-2 


1575 


18T/  +  162T^  +  53T^  +  -  36^.  +  34?^  -  32T^ 


and 


-  22Tg  -  34T,  ♦  32T/p  +  12T,,  -  8CT/  J  , 


k<6  C>  <CU  4>  +  KC°>+  C°SCS6+ 

.  1  ;  5T,  -  1ST  +  |t3-  %  +  11TX  -  ITT  +  2T7  ♦  4T 

315  L 


-  4T,  -  2T/0  -  6T//  +  8T/7 


133 


Henoe 


<A*> 


1575 


28T/  +  252T  +  80^  +  3^r 


56t  +  49T 

v5"  6 


and 


<B*> 


-  42T_  -  42T.  -  49T  +  42T  +  42T.,  -  126t,, 

/  a  t  /c>  7  7  ' jL  j 


95fr,  +  336t^  +  ie^3  -  53fr^  -  16^, 


1575 


-  5l^r6  -  56?^  -  56Tfi  -  123|t^  +  56t/£?+  56t//  ~  l68T/x] 


* 


These  expressions  can  be  further  simplified  by  setting 


(  a  *  Cyr  +  C  )  -  3(  + 


•M-  SS  V*  66  Ss  66 


(  c„  +  +  cS3)  -  3(  C,,C„  +  c„c„  +  c„c„) 


tt  22  66  33  22  33 


(  C/^+  C/y  +  “  5(  Cy^C/?  +  C/7C7?  +  C/?C7?)^ 


6Z  63  T  "/l  23  63  23 


C//  (C/4  +  "  2C^3^  +  CZz(C62+  CZ3  ~ 


/i1 


+  Cyy(C/5  +  “  2C/J  - 


A. 13  and  A. 14  then  become 


(A j  y  »  N  *>  2P  +  1  (  24a  ■+■  7b  ■+■  13c  +  d  )  } 

225  135 

>2 

(B  \  =  8P  +  8  (  4a  +  2b  +  3c  +  d  )  / 

N  W,  U5 

(k,\  o  +  _JL  (  12a  +  b  +  4o  -  2d  )  # 

6  150  90 


A.  14a 


A.  14b 


A. 15 


A. 16 
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